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Abstract 

There exist many experimental indications that final-state interactions (FSIs) may play a promi- 



> . 

, nent role not only in charmful B decays but also in charmless B ones. We examine the final-state 



rescattering effects on the hadronic B decay rates and their impact on direct CP violation. The 
On , color-suppressed neutral modes such as B^ D , K^TT^ can be substantially en- 

, hanced by long-distance rescattering effects. The direct CP- violating partial rate asymmetries in 

^ i charmless B decays to tttt/ttK and are significantly affected by final-state rescattering and their 



signs are generally different from that predicted by the short-distance approach. For example, di- 
rect CP asymmetry in B^ p^n^ is increased to around 60% due to final state rescattering effects 



D , whereas the short-distance picture gives about 1%. Evidence of direct CP violation in the decay 

' B^ K^TT^ is now established, while the combined BaBar and Belle measurements of B^ p^ir^ 



imply a 3. 60" direct CP asymmetry in the p^Tr~ mode. Our predictions for CP violation agree 
with experiment in both magnitude and sign, whereas the QCD factorization predictions (espe- 
5^ I cially for p'^Tr~) seem to have some difficulty with the data. Direct CP violation in the decay 

B~ 7r~7r'' is very small ( < 1%) in the Standard Model even after the inclusion of FSIs. Its 
measurement will provide a nice way to search for New Physics as in the Standard Model QCD 
penguins cannot contribute (except by isospin violation). Current data on ttK modes seem to 
violate the isospin sum rule relation, suggesting the presence of electroweak penguin contributions. 
We have also investigated whether a large transverse polarization in i? — > (pK* can arise from the 
final-state rescattering of D^*^Di*'^ into (j)K*. While the longitudinal polarization fraction can be 
reduced significantly from short-distance predictions due to such FSI effects, no sizable perpen- 
dicular polarization is found owing mainly to the large cancellations occurring in the processes 
B — > D*D — > (j)K* and B DgD* <j)K* and this can be understood as a consequence of 
CP and SU(3) [CPS] symmetry. To fully account for the polarization anomaly (especially the 
perpendicular polarization) observed in i? — > (pK* , FSI from other states or other mechanism, e.g. 
the penguin-induced annihilation, may have to be invoked. Our conclusion is that the small value 
of the longitudinal polarization in VV modes cannot be regarded as a clean signal for New Physics. 
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I. INTRODUCTION 



The importance of final-state interactions (FSIs) has long been recognized in hadronic charm 
decays since some resonances are known to exist at energies close to the mass of the charmed 
meson. As for hadronic B decays, the general folklore is that FSIs are expected to play only a 
minor role there as the energy release in the energetic B decay is so large that the final-state 
particles are moving fast and hence they do not have adequate time for getting involved in final- 
state rescattering. However, from the data accumulated at B factories and at CLEO, there are 
growing indications that soft final-state rescattering effects do play an essential role in B physics. 

Some possible hints at FSIs in the B sector are: 

1. There exist some decays that do not receive any factorizable contributions, for example B 
Kxoc, owing to the vanishing matrix element of the {V — A) current, (xoc|c7^(l — 75)c|0) = 0. 
Experimentally, it was reported by both Belle 01 and BaBar that this decay mode has 
a sizable branching ratio, of order (2 ~ 6) x 10~^. This implies that the nonfactorizable 
correction is important and/or the rescattering effect is sizable. Studies based on the light- 
cone sum rule approach indicate that nonfactorizable contributions to i? — > XcO^^ due to soft 
gluon exchanges is too small to accommodate the data 0, |^. In contrast, it has been shown 
that the rescattering effect from the intermediate charmed mesons is able to reproduce the 
observed large branching ratio 5]. 

2. The color-suppressed modes ^ D^*'''^tt^ have been measured by Belle 0, CLEO [3| and 
BaBar 0|. Their branching ratios are all significantly larger than theoretical expectations 
based on naive factorization. When combined with the color-allowed B D^*^ti decays, it 
indicates non- vanishing relative strong phases among various B D^*^tt decay amplitudes. 
Denoting T and C as the color-allowed tree amplitude and color-suppressed I^-emission 
amplitude, respectively, it is found that C/T ~ 0.45 exp(ibi60°) (see e.g. 'j)]), showing a 
non-trivial relative strong phase between C and T amplitudes. The large magnitude and 
phase of C/T compared to naive expectation implies the importance of long-distance FSI 
contributions to the color-suppressed internal M^-emission via final-state rescattering of the 
color-allowed tree amplitude. 

3. A model independent analysis of charmless B decay data based on the topological quark 
diagrammatic approach yields a larger value of \C/T\ and a large strong relative phase 
between the C and T amplitudes |l(J |. For example, one of the fits in leads to 
C/T = (0.46lo;3o) exp([— i(94l52)]°) which is indeed consistent with the result inferred from 
B Dtt decays. This means that FSIs in charmless B decays are as important as in S ^ Dtt 
ones. The presence of a large color-suppressed amplitude is somewhat a surprise from the 
current model calculations such as those in the QCD factorization approach and most likely 
indicates a prominent role for final-state rescattering. 



4. Both BaBar 



11( 1 and Belle |12|] have reported a sizable branching ratio of order 2 x 10 for 



the decay B^ — > n^ir^. This cannot be explained by either QCD factorization or the PQCD 
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approach and it again calls for a possible rescattering effect to induce Likewise, the 

color-suppressed decay B with the branching ratio (1.9 it 1.2) x 10 averaged from 

the Belle and BaBar measurements (see Table IV)) is substantially higher than the prediction 
based on PQCD or the short-distance approach. Clearly, FSI is one of the prominent 
candidates responsible for this huge enhancement. 

5. Direct CP violation in the decay K^n" with the magnitude Ak-k = —0.133 it 0.030 ± 

0.009 at 4.2(7 was recently announced by BaBar which agrees with the previous Belle 
measurement of jr = —0.088 it 0.035 it 0.013 at 2Aa based on a 140 fb~^ data sample 
A further evidence of this direct CP violation at 3.9a was just reported by Belle with 253 fb~^ 
of data In the calculation of QCD factorization [l^. [l^. the predicted CP asymmetry 
Aktt = (4.5]'jg'g)% is apparently in conflict with experiment. It is conceivable that FSIs via 
rescattering will modify the prediction based on the short-distance interactions. Likewise, 
a large direct CP asymmetry in the decay B^ vr+vr^ was reported by Belle Q|, but 
it has not been confirmed by BaBar The weighted average of Belle and BaBar gives 

Atttt = 0.31 ± 0.24 ji^] with the PDG scale factor of = 2.2 on the error. Again, the central 
value of the QCD factorization prediction fli has a sign opposite to the world average. 



6. Under the factorization approach, the color-suppressed mode B — > Dj^K can only proceed 



22 1 and BaBar 



via M^-exchange. Its sizable branching ratio of order 4x 10~^ observed by Belle 
[2^ will need a large final-state rescattering contribution if the short-distance H^-exchange 
effect is indeed small according to the existing model calculations. 



7. The measured longitudinal fractions for B <j)K* by both BaBar j2J] and Belle [25| are close 
to 50%. This is in sharp contrast to the general argument that factorizable amplitudes in B 
decays to light vector meson pairs give a longitudinal polarization satisfying the scaling law: 
1 — fi = 0{l/m1). ^ This law remains formally true even when nonfactorizable graphs are 
included in QCD factorization. Therefore, in order to obtain a large transverse polarization 
in S — > (pK* , this scaling law valid at short-distance interactions must be violated. The 
effect of long-distance rescattering on this scaling law should be examined. 

The presence of FSIs can have interesting impact on the direct CP violation phenomenology. 
As stressed in [3| ■, traditional discussions have centered around the absorptive part of the penguin 
graph in 6 — > s transitions [2^ and as a result causes "simple" CP violation; long-distance final 
state rescattering effects, in general, will lead to a different pattern of CP violation, namely, 
"compound" CP violation. Predictions of simple CP violation are quite distinct from that of 
compound CP violation. The sizable CP asymmetry observed in ^ K^ir^ decays is a strong 
indication for large direct CP violation driven by long-distance rescattering effects. Final state 



More recently Kagan has suggested that the penguin-induced annihilation can cause an appreciable devia- 
tion from this expectation numerically, though the scaling law is still respected formally |2(]ll27l |. However, 
it is difficult to make a reliable estimation of this deviation. 
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rescattering phases in B decays are unlikely to be small possibly causing large compound CP- 
violating partial rate asymmetries in these modes. As shown below, the sign of CP asymmetry 
can be easily flipped by long-distance rescattering effects. Hence, it is important to explore the 
compound CP violation. 

Of course, it is notoriously difficult to study FSI effects in a systematic way as it is nonperturba- 
tive in nature. Nevertheless, we can gain some control on rescattering effects by studying them in 
a phenomenological way. More specifically, FSIs can be modelled as the soft rescattering of certain 
intermediate two-body hadronic channels, e.g. B DD —>■ vrvr, so that they can be treated as the 
one-particle-exchange processes at the hadron level. That is, we shall study long-distance rescatter- 
ing effects by considering one-particle-exchange graphs. As the exchanged particle is not on-shell, 
form factors must be introduced to render the whole calculation meaningful in the perturbative 
sense. This approach has been applied to the study of FSIs in charm decays for some time |3fll. l3]]|. 
In the context of B physics, the so-called "charming penguin" contributions to charmless hadronic 
B decays have been treated as long-distance effects manifesting in the rescattering processes such 
as B ^ DgD — > Kt: 32l. 1331. l34l| . Likewise, the dynamics for final-state interactions is assumed to 
be dominated by the mixing of the final state with D^*^D^*^ in 13 5l . l3q | . Final-state rescattering 
effects were also found in B to charmonium final states, e.g. B^ K~XcO, a process prohibited in 
the naive factorization approach j^. Effects of final-state rescattering on Kit and vrvr final states 



28, 



37, 



have been discussed extensively in the literature 

The approach of modelling FSIs as soft rescattering processes of intermediate two-body states 
has been criticized on several grounds [s^. First, there are many more intermediate multi-body 
channels in B decays and systematic cancellations among them are predicted to occur in the heavy 
quark limit. This effect of cancellation will be missed if only a few intermediate states are taken 
into account. Second, the hadronic dynamics of multi-body decays is very complicated and in 
general not under theoretical control. Moreover, the number of channels and the energy release 
in B decays are large. We wish to stress that the b quark mass (~ 4.5 GeV) is not very large 
and far from infinity in reality. The aforementioned cancellation may not occur or may not be 
very effective for the finite B mass. For intermediate two-body states, we always consider those 
channels that are quark-mixing-angle most favored so that they give the dominant long-distance 
contributions. Whether there exist cancellations between two-body and multi-body channels is not 
known. Following we may assume that two-body ^ n-body rescatterings are negligible either 
justified from the 1/Nc argument or suppressed by large cancellations. We view our treatment 
of the two-body hadronic model for FSIs as a working tool. We work out the consequences of 
this tool to see if it is empirically working. If it turns out to be successful, then it will imply 
the possible dominance of intermediate two-body contributions. In other approaches such as QCD 
factorization jl^], the complicated hadronic B decays in principle can be treated systematically 
as l/mfj power corrections which vanish in the heavy quark limit. However, one should recognize 
that unless the coefficients of power corrections are known or calculable in a model independent 
manner, the short-distance picture without a control of the 1/mf, coefficients and without being 
able to include the effects of FSIs have their limitation as well. When the short-distance scenario is 
tested against experiment and failings occur such as the many examples mentioned before, then it 
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does not necessarily mean that New Physics has been discovered and an examination of rescattering 
effects can be helpful. 

The layout of the present paper is as follows. In Sec. II we give an overview of the role played by 
the final-state interactions in hadronic B decays. As a warm-up, we begin in Sec. Ill with a study of 
final-state rescattering contributions to — > Dn decays which proceed only through tree diagrams. 
We then proceed to the penguin-dominated B Ktt decays in Sec. IV, tree-dominated B ^ tttt 
decays in Sec. V and B ^ pn decays in Sec. VI. Effects of FSIs on the branching ratios and direct 
CP asymmetries are studied. Sec. VII is devoted to the polarization anomaly discovered recently 
in B ^ (j)K* decays. Conclusion and discussion are given in Sec. VIII. Appendix A gives some 
useful formula for the phase-space integration, while the theoretical input parameters employed in 
the present paper are summarized in Appendix B. 



II. FINAL-STATE INTERACTIONS 

In the diagrammatic approach, all two-body nonleptonic weak decays of heavy mesons can be 
expressed in terms of six distinct quark diagrams T, the color-allowed external W- 

emission tree diagram; C, the color-suppressed internal VF-emission diagram; £, the H^-exchange 
diagram; A, the PF-annihilation diagram; V, the penguin diagram; and V, the vertical VF-loop 
diagram. It should be stressed that these quark diagrams are classified according to the topologies 
of weak interactions with all strong interaction effects included and hence they are not Feynman 
graphs. All quark graphs used in this approach are topological with all the strong interactions 
included, i.e. gluon lines are included in all possible ways. 

As stressed above, topological graphs can provide information on final-state interactions (FSIs). 
In general, there are several different forms of FSIs: elastic scattering and inelastic scattering such 
as quark exchange, resonance formation,- • •, etc. Take the decay B^ — > D^tt^ as an illustration. 
The topological amplitudes C, £, A can receive contributions from the tree amplitude T of e.g. 
B^ D^TT^ via final-state rescattering, as illustrated in Fig. 1: Fig. 1(a) has the same topology as 
T^-exchange, while 1(b) and 1(c) mimic the internal M^-emission amplitude C. Therefore, even if the 
short-distance W^-exchange vanishes, a long-distance T^-exchange can be induced via inelastic FSIs. 



Historically, it was first pointed out in 45| that rescattering effects required by unitarity can produce 
the reaction — > IC^cj), for example, even in the absence of the Vl^-exchange diagram. Then it 
was shown in that this rescattering diagram belongs to the generic VF-exchange topology. 

Since FSIs are nonperturbative in nature, in principle it is extremely difficult to calculate their 
effects. It is customary to evaluate the lon g-d istance VF-exchange contribution. Fig. 1(a), at the 
hadron level manifested as Fig. Fig. 2(a) shows the resonant amplitude coming 

from ^ D^TT^ followed by a s-channel = O'^ particle exchange with the quark content (cn), 
which couples to D^tt^ and D^n^ . Fig. 2(b) corresponds to the t-channel contribution with a 
p particle exchange. The relative phase between T and C indicates some final-state interactions 
responsible for this. Figs. 1(b) and 1(c) show that final-state rescattering via quark exchange has 
the same topology as the color-suppressed internal VF-emission amplitude. At the hadron level. 
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(c) 

FIG. 1: Contributions to —>■ D^tt^ from the color-allowed weak decay —>■ D^Tr~ followed by 
a resonant-like rescattering (a) and quark exchange (b) and (c). While (a) has the same topology 
as the T^-exchange graph, (b) and (c) mimic the color-suppressed internal VF-emission graph. 

Figs. 1(b) and 1(c) are manifested in the rescattering processes with one particle exchange in the 
t channel jiB ] (see Fig. 3). Note that Figs. 3(a) and 2(b) are the same at the meson level even 
though at the quark level they correspond to different processes, namely, annihilation and quark 
exchange, respectively. In contrast, Fig. 3(b) is different than Fig. 3(a) in the context of the 
t-channel exchanged particle. 

For charm decays, it is expected that the long-distance M^-exchange is dominated by resonant 
FSIs as shown in Fig. 2(a). That is, the resonance formation of FSI via qq resonances is probably 
the most important one due to the fact that an abundant spectrum of resonances is known to exist 
at energies close to the mass of the charmed meson. However, a direct calculation of this diagram 
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FIG. 2: Manifestation of Fig. 1(a) as the long-distance s- and i-channel contributions to the 
VF-exchange ampHtude in — > D^tt^. The thick hne in (a) represents a resonance. 




FIG. 3: Manifestation of Figs. 1(b) and 1(c) as the long-distance t-channel contributions to the 
color-suppressed internal PF-emission amplitude in D^ir^. 



is subject to many theoretical uncertainties. For example, the coupling of the resonance to Dn 
states is unknown and the off-shell effects in the chiral loop should be properly addressed jZ^ . 



Nevertheless, as emphasized in |47l. l48l| . most of the properties of resonances follow from unitarity 
alone, without regard to the dynamical mechanism that produces the resonance. Consequently, as 
shown in 



47 



49|, the effect of resonance-induced FSIs [Fig. 2(a)] can be described in a model- 
independent manner in terms of the mass and width of the nearby resonances. It is found that the 
£ amplitude is modified by resonant FSIs by (see e.g. [I^) 

£ = e + (e^*^-- - 1) f e + ^ V (2.1) 



3 

with 

e2«'5- = l-i — , (2.2) 

m/5 — mji + ir/2 

where the reduced VT-exchange amplitude £ before resonant FSIs is denoted by e. Therefore, 
resonance-induced FSIs amount to modifying the VF-exchange amplitude and leaving the other 
quark-diagram amplitudes T and C intact. We thus see that even if the short-distance V7-exchange 
vanishes (i.e. e = 0), as commonly asserted, a long-distance Vl^-exchange contribution still can be 
induced from the tree amplitude T via FSI rescattering in resonance formation. 
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In B decays, in contrast to D decays, the resonant FSIs will be expected to be suppressed 
relative to the rescattering effect arising from quark exchange owing to the lack of the existence of 
resonances at energies close to the B meson mass. This means that one can neglect the s-channel 
contribution from Fig. 2(a). 

As stressed before, the calculation of the meson level Feynman diagrams in Fig. 2 or Fig. 
3 involves many theoretical uncertainties. If one naively calculates the diagram, one obtains an 
answer which does not make sense in the context of perturbation theory since the contributions 
become so large that perturbation theory is no longer trustworthy. For example, consider the loop 
contribution to B^ vr+vr^ via the rescattering process D^D^ vr+vr". Since B^ tt^tt^ is 
CKM suppressed relative to B^ D^D^ , the loop contribution is larger than the initial B in: 
amplitude. Because the t-channel exchanged particle is not on-shell, as we shall see later, a form- 
factor cutoff must be introduced to the vertex to render the whole calculation meaningful. 



III. B ^ Dtt decays 

The color-suppressed decays of B^ into D^*^^7r^ , D^r], D^lo and D^p^ have been observed h 
Belle and B^ decays into D^*^^it^ have been measured by CLEO 0]. Recently, BaBar 
has presented the measurements of ^ decays into D'^*^^ (tt^ , t] , m) and D^r]'. All measured color- 
suppressed decays have similar branching ratios with central values between 1.7 x 10~^ and 4.2 x 
10~^. They are all significantly larger than theoretical expectations based on naive factorization. 
For example, the measurement BiB^ D^tt^) = 2.5 x 10~^ (see below) is larger than the theoretical 
prediction, (0.58 ~ 1.13) x 10"'* j9|, by a factor of 2 ~ 4. Moreover, the three B —>■ Dtt amplitudes 
form a non-flat triangle, indicating nontrivial relative strong phases between them. In this section 
we will focus on i? — > Dtt decays and illustrate the importance of final-state rescattering effects. 

In terms of the quark-diagram topologies T, C and £, where T is the color-allowed external 
T^-emission tree amplitude, C, £ are color-suppressed internal W^-emission and T^-exchange am- 
plitudes, respectively, the B Dtt amplitudes can be expressed as 

T + 8, 

T + C, (3.1) 

^{-c + «), 

and they satisfy the isospin triangle relation 

^ D+^-) = V2A(B° ^ D\^) + A{B- ^ L»°7r-). (3.2) 

13(B°^D+7r-) = (2.76 ± 0.25) X 10"^ 

B{B- D^TT-) = (4.98 ±0.29) x 10~^ (3.3) 

= (2.5 ± 0.2) X 10"^ from the measi; 
2.9 ± 0.2 ± 0.3) X lO"'', BaBar_ 

±0.12 ±0.23) X 10" 
(2.74t°|^ ± 0.55) X 10-^ CLEO 



AiB"- 




A{B- - 


-.D\- 


a(b'- 





Using the data |5| 



and the world average B{B^ D^tt^) = (2.5 ± 0.2) x 10 ^ from the measurements 

r (2.9 ±0.2 ±0.3) X 10-^ 
B(B^ L»V°) = < (2.31 ± 0.12 ± 0.23) x 10"^, Belle M (3.4) 
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we find (only the central values for phase angles are shown here) 



C-£ 



T + £ 



Dn 



(0.41 ±0.05) e 



±i53° 



C-£ 



T + C 



Dtt 



(0.29 ± 0.03) e 



±i38° 



where we have used the relation, for example 



cos 9 



{C-£,T+£} 



2JB{B^ D+TT-) \/2B{B' 



(3.5) 



(3.6) 



to extract the strong phases. For the numerical results in Eq. (|3.5j) . we have employed the B meson 
lifetimes t{B^) = (1.536 ± 0.014) x lO'^^^ ^nd t{B+) = (1.671 ± 0.018) x lO'^^s H- 

It is known that in Cabibbo-allowed D — > PP decays ^ the topological amplitudes T,C, £ can be 
individually extracted from the data with the results 



521 



C_ 
T 



£_ 
T 



D^PP 



(0.59 ± 0.05) e 



jll4° 



(3.7) 



(0.73 ± 0.05) e-^i52° 

D^PP 

Therefore, in charm decays the color-suppressed amplitude C is not color suppressed and the W- 
exchanged amplitude is quite sizable ! The large phase of £ is suggestive of nearby resonance 
effects. The three amplitudes T, C and £ m B ^ Dtt decays can be individually determined from 
the measurements of Dtt in conjunction with the data of DfK~ and D^rf 

In the factorization approach, the short-distance factorizable amplitudes read 

.Gp 
72 



SD 



K6Kd«i(^vr)(m| - mDUF^^'iml] 



Gp 

CsD = i-j=VcbV*da2{DT:){m\ 
"^SD = i^VcbV*da2{Di:){m% 



Bit I 



[m 



(3.^ 



where ai_2 are the parameters related to the Wilson coefficients and some calculable short-distance 
nonfactorizable effects. The annihilation form factor Fq^^'^ [m^^) is expected to be suppressed 



at large momentum transfer, 



m 



Bi 



corresponding to the conventional helicity suppression. 



Based on the argument of helicity and color suppression, one may therefore neglect short-distance 
(hard) W^-exchange contributions. In the QCD factorization approach, contrary to the parameter 
oi(L>7r), a2{D'K) is not calculable owing to the presence of infrared divergence caused by the gluon 
exchange between the emitted meson and the (B^tt^) system. In other words, the nonfactorizable 
contribution to 02 is dominated by nonperturbative effects. For charmless B decays, a typical 
value of a2 ~ 0.22exp(— i30°) is obtained in the QCD factorization approach jl7| . Recall that 
the experimental value for B — > J/xl^K is \a2{J/il'K)\ = 0.26 it 0.02 |^. However, neglecting the 
V -exchange contribution, a direct fit to the Dtt data requires that 02/01 ~ (0.45 — 0.65)e^*^''° 



The question is then why the magnitude and phase of 02/01 are so different from 
the model expectation. To resolve this difficulty, we next turn to long-distance rescattering effects. 
An early effort along this direction is based on a quasi-elastic scattering model For a recent 
study of the color-suppressed B D'^*^^M decays in the perturbative QCD approach based on kx 
factorization theorem and in soft-collinear effective theory, see jsF 



and 



59| respectively. 
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D 



D 



(d) (e) (f) 

FIG. 4: Long-distance t-channel rescattering contributions to i? — > Dir. 



A. Long-distance contributions to B ^ Dtt 

Some possible leading long-distance FSI contributions to -B — > Dtt are depicted in Fig. [l] Apart 
from the Dtt intermediate state contributions as shown in Figs. 2 and 3, here we have also included 
rescattering contributions from the D*p and D*ai intermediate states. As noted in passing, the 
s-channel contribution is presumably negligible owing to the absence of nearby resonances. Hence, 
we will focus only on the t-channel long-distance contributions. For each diagram in Fig. ^ one 
should consider all the possible isospin structure and draw all the possible sub-diagrams at the 
quark level [s^]. While all the six diagrams contribute to B~ — > D^Tr~ , only the diagrams Ufa), 
lie) and iff) contribute to D+tt- and Hb) , Efd) and He) to B° D°7r°. To see this, we 

consider the contribution of Fig. 4(a) to D^n^ as an example. The corresponding diagrams 

of Fig. 4(a) at the quark level are Figs. 1(a) and 1(b). At the meson level. Fig. 4(a) contains 
Figs. 2(b) and 3(a). Owing to the wave function vr^ = {uu — dd)/y/2, Fig. 1(a) and hence Fig. 
2(b) has an isospin factor of l/\/2, while Fig. 1(b) and hence Fig. 3(a) has a factor of —1/^/2. 
Consequently, there is a cancellation between Figs. 2(b) and 3(a). Another way for understanding 
this cancellation is to note that Fig. 2(b) contributes to 8, while Fig. 3(a) to C. From Eq. 1)3. If) , 
it is clear that there is a cancellation between them. 

Given the weak Hamiltonian in the form iJw = X)i ^iQi^ where Aj is the combination of the 
quark mixing matrix elements and Qi is a T-even local operator (T: time reversal), the absorptive 
part of Fig. 4 can be obtained by using the optical theorem and time-reversal invariant weak decay 
operator Qi. From the time reversal invariance of Q (= UtQ*u\,), it follows that 

(i;out|Q|i?;in)* =^5*,(j;out|Q|i?;in), (3.9) 

3 

where Sij = (i;out|j;in) is the strong interaction S-matrix element, and we have used 
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C/rlout (in))* = |in(out)) to fix the phase convention.^ Eq. (|3.9() imphes an identity related to 
the optical theorem. Noting that 5 = 1 + iT, we find 

2Abs{i;out\Q\B;m) = ^ T/^O'; out|Q|B; in), (3.10) 

where use of the unitarity of the S'-matrix has been made. Specifically, for two-body B decays, we 
have 

AbsMipB^pm) = 2 ( nLi / ..fiL 1 (2^)' 




X <5^(r +P2 - E 9k)MipB ^ {qk})T*{piP2 ^ {qk})- (3.11) 
k=l 

Thus the optical theorem relates the absorptive part of the two-body decay amplitude to the sum 
over all possible B decay final states {qk}, followed by strong {q^} P1P2 rescattering. 

Neglecting the dispersive parts for the moment, the FSI corrections to the topological amplitudes 

are 

r = TsD, 

C = CsB + iAhs{Aa + 4b + Ac + M + ^e + Af), (3.12) 
£ = £^Y)+iAhs{Aa + + 'if). 

The color-allowed amplitude T does receive contributions from, for example, the s-channel B^ — > 
D^TT^ D^TT~ and the t-channel rescattering process B L»°7r° D+TT-. However, they 
are both suppressed: The first one is subject to 0{l/ml) suppression while the weak decay in the 
second process is color suppressed. Therefore, it is safe to neglect long-distance corrections to T. 
As a result, 

TsD + ^^SD + iAbs (4a + 4c + 4/), 

TsD + CsD + iAbs (4a + 4b + Ac + Ad + Ae + 4/), (3.13) 
-^(-C + f )sD - -^^^^ + 4d + 4e). 

Note that the isospin relation 1)3. 2(1 is still respected, as it should be. 

3 



a(b"- 




A{B- - 


-.D\- 


a(b'- 





To proceed we write down the relevant Lagrangian |6(l . I61l 

jO = -^gp^Jp+7T^^^7T-+p;,7T+^^'7^'> + pl7T-^^'7T 



2 Note that in the usual phase convention we have r|P(p)) = -\P{-p)), T\V{p,X)) = -{-)^\V{-p, X)) 
with A being the helicity of the vector meson. For the B PP, VP, VV decays followed by two-particle 
to two-particle rescatterings, the \PP) and S',£'-wave \VV) states satisfy the J7T|out (in))* = |in(out)) 
relation readily, while for \B) and P-wave \ VV) as well as \VP) states we may assign them an additional 
phase i to satisfy the above relation. We shall return to the usual phase convention once the optical 
theorem for final-state rescattering, namely, Eq. H3.10|l . is obtained. 

^ In the chiral and heavy quark limits, the effective Lagrangian of (|3.15|l can be recast compactly in terms 
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- igDDvDj D\Vn) - 2fD*Dve^uaf3{d^Vn){Dld^D*('^ - of^ d'^D^) 

+ igD'D'vDp Dyivn- + 4ifD*D'vD*l{d^V - d^V^jpl\ (3.15) 

with the convention e^i^a _ -^^ where P and are 3x3 matrices for the octet pseudoscalar and 

nonet vector mesons, respectively 



( 



P 



V 



\ K*- 



K*+\ 



Note that our phase convention on fields is fixed by 

and is different from ^6^. In the chiral and heavy quark limits, we have [g^I 



9d*d*-k 



gD*D-K 



-T-g, 9DDV = gD*D*V 



f^gv 



fo'DV 



fD*D' 



V 



^gv 



(3.16) 



(3.17) 



(3.18) 



with = 132 MeV. The parameters gy, (3 and A thus enter into the effective chiral Lagrangian 
describing the interactions of heavy mesons with low momentum light vector mesons (see e.g. ,6(1]). 
The parameter gy respects the relation gy = rapj j^^ |6fli]. We shall follow [3| to use /3 = 0.9 and 
A = 0.56 GeV^^. Instead of writing down the Feynman rules for the vertices, we work out the 
corresponding matrix elements: 

(p+(e)7r"(p2)N'C|7r°(pi)) = -igp^ne ■ {pi + P2), 
{D{p2)iT{q)\iC\D*{e,pi)) = -ign-Dn £: ■ q, 
{D*{e2,P2Mq)\iC\D*{ei,pi)) -- 



-igD'D'TTepual3£ie2 Q Pi: 



'V2 



Pgv £ • (pi + P2) 



-i2V2\gyep,^pe^;,£%pU^. 



(3.19) 



{D{p2)p{e)\iC\D{p^)) -- 

{D*{e2,P2)p{ep,q)\iC\D{pi)) -- 

Figs. Ufa) and ^b) arise from the weak decay B Dir followed by the rescattering of Dn to 
Dir. Denoting the momenta by B{pb) D{pi)t^{P2) D{p^)Ti{pi)^ it follows that the absorptive 
part of ^a) is given by ^ 



Abs{Aa) = 
of superfields 



d?pi Sp2 



{2^f5^{pB-pi-P2)A{B'^ ^ D+Tx 



{2-Kf2Ei {2'Kf2E2 



(3.14) 



where the superfield H is given by iJ = — |^(£'*^7p — 175!?), and {Vp)ba and {Ap)ba arc the matrix elements 
of vector and axial currents, respectively, constructed from Goldstone bosons. 

As noticed in the isospin factor of 1/V2 or — 1/V2 should be dropped for the intermediate state p° 



or TT 
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V2 t -mf, + iiripTp \^ m'^p 

2/ 



Pi dcos^ 1 ^ , /-ftO ^j- s F (t,mo) 



|pi|(icos6' 1 ^ ^ ^ ^^^0 F'^{t,mp 

iji + inopj. p 

where is the angle between pi and p3, k is the momentum of the exchanged p meson {k"^ = t), 
and 

t = {pi-psY = ml + ml- 2EiE3 + 2\pi\\p3\ cos 6, 

/ N (ml - ml) (mo - ml) /ooi\ 
-H^i = -(.P1-P2+P3-P4 + P1-P4+P2-P3)-— ^^V^ -■ 3.21 

m^ 

The form factor F{t, m) in Eq. (|3.2Ujl takes care of the off-shell effect of the exchanged particle, 
which is usually parametrized as 

F{t,m)= { j^2_t ) ' (3-22) 

normalized to unity at t = m^. The monopole behavior of the form factor (i.e. n = 1) is preferred 
as it is consistent with the QCD sum rule expectation [g^]. However, we shall return back to this 
point when discussing the FSI effects in — > (pK* decays (see Sec. VII. C). 
Likewise, the absorptive part of^b) is given by 

^'^(^" = \ I m m - - "''-^"^ ^ 

. p F^{t,mD*) I , kpky \ 

X igO-D-KPA — 2 [-VgD'-D7T{-P2) -gfiu + 

t-mj^, \ mjj, J 



^ |pi|dcOS^^2 n+„-A ^ (^'''"■D* 

D* 



g^^.^^A{B"^D^n-)-^^^H2, (3.23) 
_i IbirmB t — mt)* 



with 



H2 = -P2 ■ P. + (^^•P3-mi)(pi-p4-ml) ^ 



Figs. H^c)-|He) come from the rescattering process B D*(pi)p{p2) F>{p3)-K{p4^). The 
absorptive part of^c) reads 



■^"'^'^^ = II im m j: -^'^ - °"^') 



Ai,A2 



X {-i)gD*D7T ei • (-Ps) — — {-i^gpTVTT 262 • P4- (3.25) 

I — m^ 

To proceed, we note that the factorizable amplitude of i? — > V1V2 is given by 

Gp — 
A{B^V^V2) = -^VcKM{V2\{q2q3)y.A^){Vi\{qih)^_^B) 



-ifv2m2 



(el •e*2)(rnB + mi)Af^^(m^) - (e^- ) (e*2 • p J - 



+ m\ 



iepyape2'^£i PgPi 



TUB + mi 



(3.26) 
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with (qm), 



^17/^(1 - I5)q2- Therefore, 

\pi\dcos9 



Abs (4c) 



-i2 



-1 IGirmB 



t — mi 



where 



^3 



P3 - Pi 



P2 ■ P3 



{Pi ■ Pz){Pi ■ Pa) {P2 ■ P3){P2 ■ Pa) , {pi ■ P2){pi ■ P3){P2 ■ Pa) 



(3.27) 



+ 



TJln 



2 2 



{Pi ■P2){Pi -Ps) 



Pi -PA 



{Pi ■P2){P2 - Pa) 



ml 



(3.28) 



Likewise, 



Abs {Ad) 



1 [ d^pi d^p2 4 4 ^ —0 ^+ _ 

{2TT)^d^{pB-pi-p2) 2^ A{B -^D^p 



2 J (27r)32Ei {2tt)^2E2 

F^{t,mD') 



Ai,A2 



X {-i)9D*D7T£i - Pa 



t^mlT' '^^ave2-ik + P3) 



iV2 



1 1^, 



|pi|(icos 9 
_i IGnmB 



P 9v QD'D-k 



F'^it,mD) 



t — m 



D 



2A^^*(m'^) 

X fpmp (mB + mD')Af°' {ml)H4 ^-^H', 



(3.29) 



rriB + mo* 

where H4 and H'^ can be obtained from and H^, respectively, by interchanging ^3 and p4, and 



Abs (4e) 



1 f d^pi d^p2 44 ^ ^0 

(27r) 6 [pB -Pi -P2) 2^ A{B D ^p ) 



2 J (27r)32^i {2tt)'^2E2 



Ai,A2 



X igD*D*nepua(3£lP2Pl 



— 2 — i2V2Xgv epa^jjE^Psi- 



t — mf). 



-P2)'{-r + ^ 



m 



D* 



i2V2 



1 \pi\dcose F^{t,mD*) 

—TT, gv>^9D*D*TT — 2 



X fpmp {mB + mD*)Af^* {mp)H5 



t — m^* 



mB + mo* 



(3.30) 



where 



^5 
Hi 



= 2(pi •P2)(P3 - Pa) -2(pi •P3)(P2 -^4), 

= m\y{pi ■P2){P3-Pa) - {Pi ■Pa){P2 ■Pa)\ + {P2 ■ Pb){pa ■ Pb){pi -Ps) (3.31) 
+ {Pi ■ Pb){P3 ■ Pb){P2 ■ Pa) - {pi ■ Pb){P2 ■ Pb){P3 ■ Pa) - {P3 • Pb){pa ■ Pb){pi ■ P2)- 



Fig. I^f) comes from the weak decay D*^a^ fohowed by a strong scattering. We obtain 



Abs (4/) 



^ ^ "^'^"^ {2^)'5Hpb - pi - P2) E - D*+^i) (3-32) 



2 J {2-k)'^2Ei (27r)32E2 



Al,A2 



X I 



2V2gvXep,ai3£'tP3k^ 



FHt,mp) 
t — mj, + impTp 



{gGp„ + iLp^)e^J-g''- + 



m 



D* 
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where 



,3.33) 

and Y = {p2 ■ k)"^ — k'^m^- The parameters g and I appearing in Eq. p.32|l also enter into the 
strong decay amphtude of oi pir parametrized as 

A{ai ^ pn) = {gG^, + ^V^e^ (3.34) 

and hence they can determined from the measured decay rate and the ratio of D and S waves: 

N Pc |2 , "^ai"^P |,|2^ 

£^ ^ f7i {Ep-mp)g + plma^h 

S ^\Ep + 2mp)g+p^^m,,h' ^''•"'^^ 

with Pc being the cm. momentum of the /9 or tt in the ai rest frame, and 

^ ^ Pa,_Pp f ^ ^_^<^ ) . (3.36) 




Hence, Eq. (|3.32j) can be recast to 

Abs{4:f) = i2V2 — gv^fa^ma^- o , ■ -p : ^-^^6, (3.37) 

j_i IbTrmB t — m^ + impLp ms + mo* 



with 



Cg £ m?k'^\ r 

• ^ - y j YniiP-2 ■ P'i){P2 ■ Ps) + mlipi ■ p^){pi ■ Pa) + {pi ■ P2fp3 ■ Pa 

- mlml P3-Pi- {P2 ■ Pa){pi ■ Pz){Pi ■ P2) - {pi ■ P2){P2 ■ Pz){Pi ■ Pi) 

- 2g{m{p2 ■ P3 - (pi • P2){pi • Pa)) • (3.38) 

It should be emphasized that attention must be paid to the relative sign between B PP and 
B VV decay amplitudes when calculating long-distance contributions from various rescattering 
processes. For the one-body matrix element defined by (0|^7^75|P) = ifpq^, the signs of i? ^ PP 
and B VV amplitudes, respectively, are fixed as in Eqs. (|3.8() and 1)3. 26p . This can be checked 
explicitly via heavy quark symmetry (see e.g. 63]). 

The dispersive part of the rescattering amplitude can be obtained from the absorptive part via 
the dispersion relation 

Vis A{ml) = - r '^^'^^'? ds'. (3.39) 
IT Js s' — 

Unlike the absorptive part, it is known that the dispersive contribution suffers from the large 
uncertainties due to some possible subtractions and the complication from integrations. For this 
reason, we will assume the dominance of the rescattering amplitude by the absorptive part and 
ignore the dispersive part in the present work except for the decays B^ — > tt+tt" and B^ — > vr'^vr'^ 
where a dispersive contribution arising from DD irir and vrvr vrvr rescattering via annihilation 
may play an essential role. 
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B. Numerical results 



To estimate the contributions from rescattering amplitudes we need to specify various parameters 
entering into the vertices of Feynman diagrams. The on-shell strong couphng ^pTrvr is determined 
from the p ^ vrvr rate to be Qpj^T^ = 6.05 ± 0.02. The coupling gD*Dn has been extracted by CLEO 



64 1 . The parameters relevant for the D*D^*^f 



to be 17.9 lb 0.3 it 1.9 from the measured D*~^ width 
couplings are gy = 5.8, (3 = 0.9 and A = 0.56 GeV^^. 

For form factors we use the values determined from the covariant light-front model '63\. For the 
coefficients g and £ in Eq. ()3.34l) . one can use the experimental results of D/S = —0.1 it 0.028 and 
r(ai pir) = 250 — 600 MeV |50|] to fix them. Specifically, we use D/S = —0.1 and r(ai pir) = 
400 MeV to obtain g = 4.3 and ^ = 5.8 . For decay constants, we use = 132 MeV, fn = 200 
MeV, Jd* = 230 MeV and fa, = -205 MeV.^ 

Since the strong vertices are determined for physical particles and since the exchanged particle 
is not on-shell, it is necessary to introduce the form factor F{t) to account for the off-shell effect 
of the t-channel exchanged particle. For the cutoff A in the form factor F{t) [see Eq. ()3.22|l ]. A 
should be not far from the physical mass of the exchanged particle. To be specific, we write 

A = mexc + r/AQCD, (3.40) 

where the parameter ij is expected to be of order unity and it depends not only on the exchanged 
particle but also on the external particles involved in the strong-interaction vertex. Since we do 
not have first-principles calculations of rj, we will determine it from the measured branching ratios 
given in ()3.3|) and 1)3. 4|) . Taking Aqcd = 220 MeV, we find r] = 2.2 for the exchanged particles D* 
and D and r] = 1.1 for p and vr. As noted in passing, the loop corrections will be larger than the 
initial B Dn amplitudes if the off-shell effect is not considered. Although the strong couplings are 
large in the magnitude, the rescattering amplitude is suppressed by a factor of F'^{t) ~ m^AqQ^/t^. 
Consequently, the off-shell effect will render the perturbative calculation meaningful. 
Numerically, we obtain (in units of GeV) 

A(B° L»+7r") = i6.41 X 10"^ - (0.90 + iOm) x lO"'^, 
A(B° DV) = -i0.91 X lO"'^ - 1.56 x 10""^, 

A{B- D^TT') = i7m X lO"'^ + (1.30 - i0.09) x lO"'^ (3.41) 

for ai = 0.90 and 02 = 0.25,^ where the first term on the right hand side of each amplitude 
arises from the short-distance factorizable contribution and the second term comes from absorptive 
parts of the final-state rescattering processes. Note that the absorptive part of D^tt^ and 

B^ D^TT^ decay amplitudes are complex owing to a non-vanishing p width. To calculate the 



^ The sign of fa^ is opposite to /^r as noticed in [ssf. 

^ As mentioned before, a2{DTr) is a priori not calculable in the QCD factorization approach. Therefore, 
we choose 02 (-Ptt) w 0.25 for the purpose of illustration. Nevertheless, a rough estimate of a2{D'K) has 
been made in |f 7l | by treating the charmed meson as a light meson while keeping its highly asymmetric 
distribution amplitude. It yields 02(071) « 0.25 exp(— i40°). 
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branching ratios, we will take into account the uncertainty from the Aqcd scale. Recall that 
Aqcd = 217t|f MeV is quoted in PDG 50]. Therefore, allowing 15% error in Aqcd, the flavor- 
averaged branching ratios are found to be 



B(B° - 




= (3.lt°:0) X 10-3 


(3.2 X 10" 




B{B- - 


-.D\-) 


= (5.ol°:J) X 10-3 


(4.9 X 10" 




13(B° - 




= {2.5tl:l) X 10-4 


(0.6 X 10" 





(3.42) 

where the branching ratios for decays without FSIs are shown in parentheses. We see that the D^ir^ 
mode is sensitive to the cutoff scales Ad* and A^i, whereas the other two color-allowed channels are 
not. It is clear that the naively predicted B{B^ D^n^) before taking into account rescattering 
effects is too small compared to experiment. In contrast, D^tt^ and B^ D^tt^ are almost 

not affected by final-state rescattering. Moreover, 

'0.28e-^48°, £ 



C_ 

r 



Dtt 



0.20, T 



f0.14e^96°, withFSI 
= < (3.43) 
Dn I 0, without FSI 



where we have assumed negligible short-distance VF-exchange and taken the central values of the 
cutoffs. Evidently, even if the short-distance weak annihilation vanishes, a long-distance contribu- 
tion to VF-exchange can be induced from rescattering. Consequently, 



C-£ 



T + £ 



0.40 e-*^^, C-£ 



Dtt 




(3.44) 

Dn I 0.17, without FSI 



0.20, T + C 

Therefore, our results are consistent with experiment (|3.5|) . This indicates that rescattering will 
enhance the ratios and provide the desired strong phases. That is, the enhancement of i? = 
(C — £)/{T + £) relative to the naive expectation is ascribed to the enhancement of color-suppressed 
VF-emission and the long-distance M^-exchange, while the phase of R is accounted for by the strong 
phases of C and £ topologies. 

Because the rescattering long-distance amplitudes have the same weak phases as the short- 
distance factorizable ones, it is clear that there is no direct CP violation induced from rescattering 
FSIs. 

Since the decay 5^ — > DfK^ proceeds only through the topological VF-exchange diagram, the 
above determination of £ allows us to estimate its decay rate. From Eq. (|3.43j) it follows that 
\£/{T + )P = 0.020 in the presence of FSIs. Therefore, we obtain 

-^-n = 0.019 , (3.45) 

r(5° ^ D+TT-) 

which is indeed consistent with the experimental value of 0.014 it 0.005 50']. The agreement will 
be further improved after taking into account SU(3) breaking in the £ amplitude of DfK^. 



IV. B^nK DECAYS 

The penguin dominated B irK decay amplitudes have the general expressions 
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TABLE I: Experimental data for CP averaged branching ratios (top, in units of 10 ^) and CP asym- 
metries (bottom) for B ttK j5C 





Mode 


BaBar 


Belle 


CLEO 


Average 


B- 


— > vr ii 


26.0 ± 1.3 ± 1.0 


22.0 ± 1.9 ± 1.1 


18 8+3.7+2.1 


24.1 ± 1.3 




^ 71+ K- 


17.9 ±0.9 ±0.7 


18.5 ± 1.0 ±0.7 


IO r,+2.3+1.2 

i».U_2,i_o.9 


18.2 ±0.8 


B 




12.0 ±0.7 ±0.6 


12.0±1.3tJ| 


1 9 q+2.4+1.2 
i^.y_2. 2-1.1 


12.1 ±0.8 


B 




11.4 ±0.9 ±0.6 


11.7±2.3tJ:3 


-, 9 Q+4.0+1.7 
J-^-0-3. 3-1.4 


11.5 ± 1.0 






-0.087 ±0.046 ±0.010 


0.05 ±0.05 ±0.01 


0.18 ±0.24 ±0.02 


-0.020 ± 0.034 






-0.133 ±0.030 ±0.009 


-0.101 ±0.025 ±0.005 


-0.04 ±0.16 ±0.02 


-0.113 ±0.019 






0.06 ±0.06 ±0.01 


0.04 ±0.05 ±0.02 


-0.29 ±0.23 ±0.02 


0.04 ±0.04 




A-ttO Ks 


-0.06 ±0.18 ±0.06 


-0.12 ±0.20 ±0.07 




-0.09 ±0.14 






0.35t{]|^ ± 0.04 


0.30 ±0.59 ±0.11 




'^4+0-27 

'-'••5^-0.29 



A(B^ ^ TT+K-) = T + V + ^P^w + ^A, 

^T,^-^) = -L(C-P±Pew±^7'ew-^a), 
A{B' ^ TT-lf) =V- ^P^w +A + Va, 



(4.1) 



A{B- 



_L(T + C ± P ± Pew ± Ir^w + ^ + T^a), 



where Vew and Pgw color- allowed and color-suppressed electroweak penguin amplitudes, re- 
spectively, and Va is the penguin-induced weak annihilation amplitude. The decay amplitudes 
satisfy the isospin relation 



A{b" TT+K~) + V2A{b" TT^W) = -A{B- tt~k") + V2A{B- n^R- 
Likewise, a similar isospin relation holds for charge conjugate fields 

A{B° ^ TT-K+) + V2A{B'^ ^ 7r°:^°) = -A{B+ ^ tt+K^) ± V2A{B+ ^ tt^R- 



(4.2) 



(4.3) 



In the factorization approach |66l . l67t | , 

TsD = i^iKai, CsD = K2A„a2, VsD = Xu{al ± a'^r^) + Xc{al ± agr^) 



o 3 
^EW = ^f^2iXu + Ac)(-a7 ± ag), V'i^ = -Ki \u{a^r^ + <o) ± Xda^r^ + a^o) 



with Ag = V*^Vgb, = 2mj^/[mbifi){ms ± mg)(/i)], and 

= i9Zf^Ft{rnl){ml - ml), K2 = ^^UFi^ {ml){ml - ml). 



(4.4) 



(4.5) 



The parameters a^'^ can be calculated in the QCD factorization approach [17|. They are basically 
the Wilson coefficients in conjunction with short-distance nonfactorizable corrections such as vertex 
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(a) 



(b) 




(c) 




(d) 



FIG. 5: Long-distance t-channel rescattering contributions to S — > Kir. 



corrections and hard spectator interactions. Formally, ai{i 6,8) and uq^s^x should be renormal- 
ization scale and scheme independent. In practice, there exists some residual scale dependence in 
ai{fj,) to finite order. At the scale fi = 2.1 GeV, the numerical results are 



ai = 0.9921 + i0.0369, 02 = 0.1933 - i0.1130, 



03 



-0.0017 + ^0.0037, 



a 



u 
4 

u 



-0.0298 - i0.0205, 
-0.0586 -«0.0188, 
= (45.0-^5.2) x 10"^ 



^6 



a 



u 
10 



(-58.3 + i86.1) X 10" 



-0.0375 - i0.0079, 
-0.0630 - i0.0056, 
4 = (44.2 + i3.1) X 10"^ ag 
= (-60.3 + ^88.8) x 10"^ 



05 = 0.0054 - i0.0050, 

5 



a-i = «5.4 X 10 



-(953.9 + ^24.5) x 10"^ 

(4.6) 



For current quark masses, we use mb{mb) = 4.4 GeV, mc{mb) = 1.3 GeV, ms(2.1GeV) = 90 MeV 
and mq/nis = 0.044. 



coefficients a^'^ leads to 


B{B- 


IT K°)sD = 


17.8 X 10" 




-^7r+K~)sB = 


13.9 X 10~ 


B{B' 




9.7 X 10"^ 


B(B° 


- ^'k\b = 


6.3 X 10"^ 


\Vcb\ = 


0.041, \Vub/Vcb\ 


= 0.09 and 



/iSD 

aSB _ 



= 0.01, 
= 0.04, 

0.08, 

-0.04. 



(4.7) 



From Table it is clear that the predicted branching ratios are slightly smaller than experiment 
especially for tt^IC^ where the prediction is about two times smaller than the data. Moreover, 

the predicted CP asymmetry for ir^K^ is opposite to the experimental measurement in sign. By 
now, direct CP violation in B^ tt^K^ has been established by both BaBar and Belle [l5| . 

i,4>2t (f's) with a = (/)2, /9 = 01, 7 = 03 are in common usage for the angles 



^ The notations (a,/?, 7) and 
of the unitarity triangle. 
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Leading long-distance rescattering contributions to B 
amplitudes are 

'■^ \pi\dcos0 



Ktt are shown in Fig. [51 The absorptive 



Abs (5a) 
Abs (56) 
Abs (5c) 



1 IGirmB 
^ \pi\dcos'^ 

-1 



gK'Kn gK*Kri'A{B^ 



m 



K* 



D+DZ 



J 2 



m 



D* 



Pi|dcos 9 



■ gD*Dn gO'DK ' 



F\t,mD) 



m 



D 



{tub + mz). )^f ^* (ml,*) Hi 



2Af^*( 



rriB + mo* 



-Hi 



Abs (bd) 



Gf 

i—^VcbVcsfDimD 



(ms + mz).)^f^* iml),)J5 



^ |^Ji|dcos6' F'^{t,mD* 

gD*D*7r gD*DK 



t 



m 



D* 



2ABD' 



{ml,] 



TUB + TTlD* 



with 



J2 



-P3 ■ P4 



{pi ■ P3 - ml){p2 ■ Pa - ml) 



m 



(4.9) 



D* 



and J5, J5 can be obtained from and H'^, respectively, under the replacement of ps 
Hence, 

> 7r"^°)sD + iAbs{5a + 56 + 5c + 5d), 
'Jt'^K~)sb + iAbs{5a + 56 + 5c + 5d), 



Pa- 



A{B- - 




= a{b- 


A{^- 


* vr+ET") 


= a(b' 


A{B' - 




= a{b- 


a{b'- 


ttOk") 


= a{b' 



7rOi^-)sD + 



V2 



Abs{ba + bb + bc + bd), 



^^6s(5a + 56 + 5c + 5d). 
V2 



(4.10) 



As noticed before, the rescattering diagrams 5(b)-5(d) with charm intermediate states contribute 
only to the topological penguin graph. Indeed, it has been long advocated that charming-penguin 
long-distance contributions increase significantly the B — > Ktt rates and yield better agreement 
with experiment For a recent work along this line, see [3^ . 

To proceed with the numerical calculation, we use the experimental value A{B^ K^v') — 
(— 7.8-|-i4.5) X 10^^ GeV with the phase determined from the topological approach For strong 
couplings, gK*K-K defined for K*^ —>■ K^tt^ is fixed to be 4.6 from the measured K* width, and 
the K*Kr]' coupling is fixed to be 



gK*+-yK+ri' 



1 1 

^gK*+->K+7T0 — —^gK* 



(4.11) 



by approximating the t]' wave function by 7]' = -^[uu + dd + 2ss). A fit to the measured decay 
rates yields = 0.69 for the exchanged particle D or D* , while rjx* is less constrained as the 

exchanged K* particle contribution [Fig. 5(a)] is suppressed relative to the D or D* contribution. 
Since the decay B DDg is quark-mixing favored over B — > Dtt, the cutoff scale here thus should 
be lower than that in i? — > Dt: decays where = 2.2 as it should be. As before, we allow 15% 
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TABLE II: Comparison of experimental data and fitted outputs for CP averaged branching ratios 
(top, in units of 10^^) and the theoretical predictions of CP asymmetries (bottom) for B ttK. 

Mode Expt. SD SD+LD 

B{B- tt-k") 24.1 ± 1.3 17^8 23.3^^;^ 

^^+T^-\ ieoj-nc IQO inQ+5-0 



B{W^TT+K-) 18.2 ±0.8 13.9 19.3_3 ^ 

B{B-^7T^K~) 12.1 ±0.8 9.7 12.5lfj 

B(B° -n^Tf) 11.5 ± 1.0 6.3 9.ll?| 

A^-Ks -0.02 ±0.03 0.01 0.024i[j:[Joi 

A^+K- -0.11 ±0.02 0.04 -0.14lg:[]^ 

A^OK- 0.04 ±0.04 0.08 -O.llig;!]^ 

A^OKs -0.09 ±0.14 -0.04 0.03ltJ]:[!?^ 

S^OKs 0.34 ±0.28 0.79 0.78 ±0.01 



uncertainty for the Aqcd scale and obtain the branching ratios and CP asymmetries as shown 
in Table ^ where use has been made of fo* ~ fOs = 230 MeV and the heavy quark symmetry 
relations 



9D*DsK = \ -9D*DTr, 9DtDK = \ -gD*D-K- (4-12) 

V rn,£) V "^D* 

It is found that the absorptive contribution from the Kr]' intermediate states is suppressed relative 
to D^*^D^*'^ as the latter are quark-mixing-angle most favored, i.e. B D^*^D^*^ ^ B ^ Kr]' . 
Evidently, all ttK modes are sensitive to the cutoffs Ad* and Ad- We see that the sign of the direct 
CP partial rate asymmetry for B^ tt^K^ is flipped after the inclusion of rescattering and the 
predicted A^^+x- = — 0.14;;'Iq q3 is now in good agreement with the world average of —0.11 ± 0.02. 
Note that the branching ratios for all ttK modes are enhanced by (30 ~ 40)% via final-state 
rescattering. 

Note that there is a ™le relation 

2Ar(7r°i^-) - Ar(7r-K°) - AT{tt+K-) + 2AT{tt°K°) = 0, (4.13) 

based on isospin symmetry. Hence, a violation of the above relation would provide an important 
test for the presence of electroweak penguin contributions. It is interesting to check how good the 
above relation works in light of the present data. We rewrite it as 

A^+k-B{tt+K-) - 2 A ,j^B{TT^lf) r-j-o 
2A^OK-B{TT'^K-)-A^_-^oBi7r-lf) ^B-' 

The left hand side of the above equation yields 0.07 ± 2.28, while the right hand side is 0.923 ± 
0.014 0. Although the central values of the data seem to imply the need of non-vanishing 
electroweak penguin contributions and/or some New Physics, it is not conclusive yet with present 
data. For further implications of this sum-rule relation, see 
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TABLE III: Experimental data for CP averaged branching ratios (top, in units of 10 ^) and 



CP asymmetries (bottom) for B tttt 



Mode 


BaBar 


Belle 


CLEO 


Average 


B 7T^7T~ 


4.7 ±0.6 ±0.2 
1.17 ±0.32 ±0.10 


4.4 ±0.6 ±0.3 

00+0.44+0.22 
^■"-•^-0.48-0.18 


A c+l. 4+0.5 
^•'^-1.2-0.4 

< 4.4 


4.6 ±0.4 
1.51 ±0.28 


B- 7r-7r° 


5.8 ±0.6 ±0.4 


5.0 ± 1.2 ±0.5 


A f.+l. 8+0.6 
^•"-1.6-0.7 


5.5 ±0.6 


A + 


0.09 ±0.15 ±0.04 


0.58 ±0.15 ±0.07 




0.31 ±0.24" 




-0.30 ±0.17 ±0.03 


-1.00 ±0.21 ±0.07 




-0.56 ±0.34'^ 




0.12 ±0.56 ±0.06 


0.43 ± 0.5ll|]:J^ 




0.28 ±0.39 


At^-t^O 


-0.01 ±0.10 ±0.02 


-0.02 ±0.10 ±0.01 




-0.02 ±0.07 



"Error-bars of ^^+^- and S',r+ir- are scaled by = 2.2 and 2.5, respectively. When taking into account 
the measured correlation between A and S, the averages are cited by Heavy Flavor Averaging Group to 
be ^TT+TT- = 0.37 ± 0.11 and S^+t^- — —0.61 ± 0.14 with errors not being scaled up. 



V. B ^TTir DECAYS 



A. Short-distance contributions and Experimental Status 



The experimental results of CP averaged branching ratios and CP asymmetries for charmless 
B TTTT decays are summarized in Table Hill For a neutral B meson decay into a, CP eigenstate 
f, CP asymmetry is defined by 

T{B{t)^f) + TiBit)^f) f ^ ' f ^ ^ J 

where Am is the mass difference of the two neutral B eigenstates, Sf is referred to as mixing-induced 
CP asymmetry and Af = —Cf is the direct CP asymmetry. Note that the Belle measurement jl2| 
gives an evidence oi CP violation in B^ vr^vr" decays at the level of 5.2 standard deviations, 
while this has not been confirmed by BaBar 120]. As a result, we follow Particle Data Group 
to use a scale factor of 2.2 and 2.5, respectively, for the error-bars in A^^+t^- and 5'^+7r- (see Table 
1111(1 . The CP-violating parameters Cf and Sf can be expressed as 



i + |; 

where 



q A{b' ^ /) 



For vrvr modes, q/p = e'"'^^ with sin 2/3 = 0.726 ± 0.037 |66 
The B TTTT decay amplitudes have the general expressions 

A(B° ^ vr+TT-) = T + V + ^P^w + £ + rA + V, 
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A{B- 



vr TT 



-^(C-P + Pew + ^^ew 



E-Va-V), 



(5.4) 



1 



_(T + c + Pew + 7^ew), 



where V is the vertical VF-loop topological diagram (see Sec. II). In the factorization approach, the 



66 



quark diagram amplitudes are given by 



67| 



-pSD 
' EW 



3 
2 



K(A„ + Ac)(-a7 + a9), Pew 



3 



A„« + <r;j) + X.ial + agr-) 
A„(a^;j + a5'o) + Ae(agr;j + a5o; 



(5.5) 



where A„ 



and the chiral enhancement factor 



2mj^/[mb{fi){mu + md){f^)] arises from the (5 — P){S + P) part of the penguin operator Qq. The 
explicit expressions of the weak annihilation amplitude £ in QCD factorization can be found in 



exp 



It is easily seen that the decay amplitudes satisfy the isospin relation 



AiB" 



IT ' TT 



V2 



A{B- TT-TT^)+A{B 



Note that there is another isospin relation for charge conjugate fields 



TT TT 



V2 



A{B^ 



TT TT 



(5.6) 



(5.7) 



Substituting the coefficients a"'^ given in Eq. (|4.6j) into Eqs. (|5.4j) and (|5.5j) leads to^ 



- 


> vr+vr )sD 


= 7.6 X 10"^ 




= -0.05, 




B(B° - 


vr TT jsD 


= 2.7 X 10"^ 


.SD 


= 0.61, 


(5 


B{B- - 


TT TT )SD 


= 5.1 X 10"^ 


4SD 


= 5 X 10"^ 





To obtain the above branching ratios and CP asymmetries we have applied the form factor 



Fq{0) = 0.25 as obtained from the covariant light-front approach j63i | and ne glec ted the W- 



annihilation contribution. Nevertheless, the numerical results are similar to that in [1^ (see Tables 
9 and 10 of 0]) except for .4^-^o for which we obtain a sign opposite to 

We see that the predicted vr+vr" rate is too large, whereas vr'^vr'' is too small compared to 
experiment (see Table inTjl . Furthermore, the predicted value of direct CP asymmetry, A^^^_ = 
—0.05, is in disagreement with the experimental average of 0.31 it 0.24. ^ In the next subsection, 
we will study the long-distance rescattering effect to see its impact on CP violation. 

In the topological analysis in 10], it is found that in order to describe B vrvr and B ttK 
branching ratios and CP asymmetries, one must introduce a large value of \CJT\ and a non-trivial 
phase between C and T. Two different fit procedures have been adopted in 10| to fit vrvr and ttK 



^ Charge conjugate modes are implicitly included in the flavor-averaged branching ratios throughout this 
paper. 

^ In the so-called perturbative QCD approach for nonlcptonic B decays ||69j, the signs of tt+tt^ and A'^7r+ 
CP asymmetries are correctly predicted and the calculated magnitudes are compatible with experiment. 
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(a) 



(b) 





(c) 



(d),(e) 



FIG. 6: Long-distance t-channel rescattering contributions to i? — > tttt. Graphs (d) and (e) corre- 
spond to the exchanged particles D and D* , respectively. 





(a) 



(b) 



FIG. 7: Contributions to i? —>■ tt^tt^ from the color-allowed weak decay B —>■ vr+vr^ followed 
by quark annihilation processes (a) and (b). They have the same topologies as the penguin and 
VF-exchange graphs, respectively. 



data points by assuming negligible weak annihilation contributions. One of the fits yields^'' 

C 



7 = i65tiir, 



T 



= (0.46t°:«)exp[-i(94t^3)]o. (5.9) 

The result for the ratio C/T in charmless B decays is thus consistent with ()3.5|1 extracted from 
B Dtt. 



B. Long-distance contributions to i? ^ vrvr 

Some leading rescattering contributions to B ^ tttt are shown in Figs. inta)-|nte). While all 
the five diagrams contribute to B^ vr+vr" and 5^ — > vr^vr'^, only the diagrams a) and Efb) 
contribute to B~ iT~7r^. Since the it~tt^ final state in B~ decay must be in / = 2, while the 



The other fit yields \C/T\ — 1.43to3i which is unreasonably too large. This ratio is substantially reduced 
once the up and top quark mediated penguins are taken into account \10^ . 
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intermediate DD state has / = or / = 1, it is clear that B 



TT TT cannot receive rescattering 



from B~ — > 7r~7r''. It should be stressed that all the graphs in Fig. El contribute to the 

penguin amplitude. To see this, let us consider Fig. [7| which is one of the manifestations of Fig. 
ini^a) at the quark level. (Rescattering diagrams with quark exchange are not shown in Fig. [7|) 
It is easily seen that Figs. [Tfa) and[7fb) can be redrawn as the topologies V and £, respectively. 
Likewise, Figs, intc)-(e) correspond to the topological V. Consequently, 



Hence, 



11 



A{B- 



TT TT 



vr TT 



T = TsD, 

C = CsD + (6a + 66), 

£ = SsB + i-^bs {6a + 6b), 

V = 'PsTi + iAhs{6a + 6h + 6c + 6d + 6e). 

= vr+7r")sD + 2i^6s (6a + 66) + i^6s(6c + 6d + 6e), 



(5.10) 



A{W 
A{B- 



1 



'^"'''""jsD H — ■j=iAbs (6a + 66 + 6c + 6d + 6e), 
v2 

vr~vr°)sD + -^iAhs (6a + 66), 



(5.11) 



The absorptive parts of B 
states in Fig. |H1 are 

^6s(6a) = 
^6s(66) = 



^h{Pi)^h{P2) ^(P3)'^(j'4) with Ml 2 being the intermediate 



1 \pi\dcOSe 2 



1 IGTrniB 

'7f 



vr vr 



t — + inipTp 



{tub + mp)Af^{mmHs - ^^^^AH'. 



Abs (6c) 
Abs {6d) 

Abs (6e) 



1 \pi\dcos9 2 

9D*DnA{B 



1 IGirmB 
. Gp 



rriB + rrip 



m 



(5.12) 



D* 



^ \pi\dcos9 2 F'^{t,n^D) 
-1 16iTmB 



dD'Dn 



m 



D 



2A. 



BD* 



BD*t^2 \zj ^^2 {'^'d*> j^i 



Of 

V2 J-i Wttitib 



rriB + niD* 
1 \pi\dcos6 2 F'^{t,mD* 



■ 9d*D*tt 



t 



m 



D* 



{mB + mD*)Af^'{m%,)J5 ^ 



(ml.) ^, 



rriB + mo* 



-Jk 



In a similar context, it has been argued in ^31] that the analogous Figs. Ela)-(b) do not contribute to the 
decay — > tt'^tt*' owing to the cancellation between the quark annihilation and quark exchange diagrams 
at the quark level. However, a careful examination indicates that the quark annihilation processes gain 
an additional factor of 2 as noticed in [s^- Consequently, the intermediate tt+tt" state does contribute to 
-B° 7r°7r° and D'^ 7r°7r° via final-state rescattering. 
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Since in SU(3) flavor limit the final-state rescattering parts of Figs. intc)-(e) should be the same 
as that of Figs, intb)-(d), it means that the cutoff scales A^) and A^i* appearing in the rescattering 
diagrams of B ^ vrvr are identical to that in S — > ttK decays in the limit of SU(3) symmetry. 
Allowing 20% SU(3) breaking and recalling that r?£)(*) = 0.69 in 5 — > ttK decays, we choose 
VdM — = 0-83 for S ^ vrvr decays. Because the decay B^ D^D^ is Cabibbo suppressed 
relative to B^ D^D^ , it is evident that the absorptive part of the charm loop diagrams, namely. 
Figs. IHc)-(e), cannot enhance the vr^vrO rate sizably. However, since the branching ratio of the 
short-distance induced B^ — > vr+vr", namely 7.6 x 10^^ [see Eq. 1)5. 8() ]. already exceeds the measured 
value substantially, this means that a dispersive part of the long-distance rescattering amplitude 
must be taken into account. This dispersive part also provides the main contribution to the 
rate. There is a subtle point here: If the dispersive contribution is fixed by fitting to the measured 
B — > vrvr rates, the value of order 0.2 is excellent and the predicted direct CP violation for vr"'"vr~ 
agrees with experiment. Unfortunately, the calculated mixing- induced CP- violating parameter S 
will have a wrong sign. Therefore, we need to accommodate the data of branching ratios and the 
parameter S simultaneously. We find (in units of GeV) 

Vis A = 1.5 X IQ-^VcbV^a - 6.7 x 10" V^bT/*^, (5.13) 

for r^TT = ?7p = 1.2. It should be stressed that these dispersive contributions cannot arise from the 
rescattering processes in Figs. inta)-(e) because they will also contribute to B ^ Kit decays via 
SU(3) symmetry and modify all the previous predictions very significantly. As pointed out in |38l |. 
there exist vrvr — > vrvr and DD — > vrvr meson annihilation diagrams in which the two initial quark 
pairs in the zero isospin configuration are destroyed and then created. Such annihilation diagrams 
have the same topology as the vertical VF-loop diagram V as mentioned in Sec. II. Hence, this FSI 
mechanism occurs in vr^vr" and but not in the i^vr modes. That is, the dispersive term given 
in Eq. 1)5.13(1 does not contribute to i? ^ Ktt decays. 

The resultant amplitudes are then given by (in units of GeV)^^ 

(2.03 + i2.02) X 10"^ + (1.40 - iiO.77) x 10"^ - ilM x 10~^ 
(-2.21 + i2.04) X 10"^ - (2.21 - i2.05) x 10"^ - il.64 x 10"^ 
(-6.15 + i3.46) X 10"^ + (7.89 - i2.73) x 10"^ - «1.07 x 10"^ 
(3.32 + il.ll) X 10"^ + (3.32 + il.ll) x 10"^ - iil.07 x 10"^ 
(2.05 + ilM) X 10"^ + (2.07 - i2.73) x 10"^ 

(-1.90 + il.34) X 10"^ - (1.90 - il.34) x 10"^ (5.14) 

where the numbers in the first parentheses on the right hand side are due to short-distance con- 
tributions, while the second term in parentheses and the third term arise from the absorptive 
and dispersive parts, respectively, of long-distance rescattering. (The dispersive contribution to 
B~ — > 7r~7r'^ is too small and can be neglected.) Form the above equation, it is clear that the 



Note that B PP and B VV decay amplitudes are accompanied by a factor of i in our phase 
convention [see Eq. ()3.8|l]. 



AiB^- 


> vr^vr ) = 


A{B^ - 


> vr"''vr~) = 


a(b'- 


-^vrV) = 


A{B° - 


-^^V) = 


A{B- - 


-^vr-vr°) = 


A{B+ - 


-.vr+vrO) = 



26 



TABLE IV: Same as Table HII except for B vrvr decays. 



Mode Expt. SD SD+LD 

13(B^ vr+vr-) 4.6 ± 0.4 7^6 le^? 

^(S°^7r0^0) 1.5 ±0.3 0.3 l.bto.l 

B{B-^TT-TT^) 5.5 ±0.6 5.1 5.4 ±0.0 

A^+^- 0.31 ±0.24 -0.05 0.35l[j:i^ 

5^+^- -0.56 ±0.34 -0.66 -0.16lg:J^ 

^^0^0 0.28 ±0.39 0.56 -0.30l[j:|54 

A^-^o -0.02 ± 0.07 5 X 10-5 -0.009l[]:g[}? 



long-distance dispersive amplitude gives the dominant contribution to vr'^vr'' and contributes de- 
structively with the short-distance vr+vr- amplitude. The corresponding CP averaged branching 
ratios and direct CP asymmetries are shown in Table ITVl where the ranges indicate the uncertain- 
ties of the QCD scale Aqcd by 15%. We only show the ranges due to the cut-off as numerical 
predictions depend sensitively on it. It is evident that the ■K~Tr^ rate is not significantly affected by 
FSIs, whereas the tt+tt- and vr'^vr'' modes receive sizable long-distance corrections.^^ 

Since the branching ratio of — > D^D~ is about 50 times larger than that of B^ it^tt~ , 
it has been proposed in [s^ that a small mixing of the vrvr and DD channels can account for the 
puzzling observation of B — > However, the aforementioned dispersive part of final-state 

rescattering was not considered in and consequently the predicted branching ratio B{B^ 
vr+vr-) ~ 11 X 10~^ is too large compared to experiment. In the present work, we have shown 
that it is the dispersive part of long-distance rescattering from DD vrvr that accounts for the 
suppression of the tt+tt- mode and the enhancement of tt^tt^. 

As far as CP rate asymmetries are concerned, the long-distance effect will generate a sizable 
direct CP violation for ■7T~^7r~ with a correct sign. Indeed, the signs of direct CP asymmetries in 
B vrvr decays are all flipped by the final-state rescattering effects. It appears that FSIs lead 
to CP asymmetries for vr^vr^ and tt^tt^ opposite in sign. The same conclusion is also reached 
in the vrvr and DD mixing model Issl and in an analysis of charmless B decay data [s^ based 
on a quasi-elastic scattering model (see also 0]). However, the perturbative QCD approach 
based on the factorization theorem predicts the same sign for tt^vt- and CP asymmetries, 
namely, A^^+t^- = 0.23±0.07 and At^Qt^o = 0.30±0.10 |6^. Moreover, a global analysis of charmless 
B —I- PP decays based on the topological approach yields A^^+t^- ~ 0.33 and .A^o^o ~ 0.53 llo| (see 



A different mechanism for understanding the tttt data is advocated in .3_8J where FSI is studied based on 
a simple two parameter model by considering a quasi-elastic scattering picture. A simultaneous fit to the 
measured rates and CP asymmetries indicates that tt+tt^ is suppressed by the dispersive part of tttt — s- tttt 
rescattering, while tt^tt" is enhanced via the absorptive part. This requires a large SU(3) rescattering 
phase difference of order 90° ~ 100°. The fitted tt-tt" rate is somewhat below experiment. In the present 
work, it is the dispersive part of the rescattering from DD to tttt via annihilation that plays the role for 
the enhancement of tt'^tt'^ and suppression of Tr+Tr". 



27 



also [t^] for the same conclusion reached in a different context). It is worth mentioning that the 
isospin analyses of i? — > tttt decay data in 7l| and 72] all indicate that given the allowed region 
of the vr^vr'^ branching ratio, there is practically no constraints on vr^^Tr" direct CP violation from 
BaBar and/or Belle data. And hence the sign of AtjO^^q is not yet fixed. At any rate, even a sign 
measurement of direct CP asymmetry in vr'^vr'' will provide a nice testing ground for discriminating 
between the FSI and PQCD approaches for CP violation. 

It is interesting to notice that based on SU(3) symmetry, there are some model-independent 
relations for direct CP asymmetries between vrvr and ttK modes [t^I 

Ar(^+^-) = -Ar(7r+i^-), Ar(^°7r°) = -Ar(7r°K°), (5.15) 

with Ar(/i/2) = T{B /1/2) — r(i? — > /1/2). From the measured rates, it follows that 

It appears that the first relation is fairly satisfied by the world averages of CP asymmetries. Since 
direct CP violation in the Tr^K~ mode is now established by B factories, the above relation implies 
that the vr+vr^ channel should have A^^+t^- ~ 0(0.4) with a positive sign. For the neutral modes 
more data are clearly needed as the current experimental results have very large errors. 
Finally we present the results for the ratios C/T, £/T and V/T: 



C_ 
T 



0.36 e 



-i55° 



£_ 

r 



0.19e 



-i85° 



V 

r 



0.56 e 



i230° 



(5.17) 



to be compared with the short-distance result C/T\^^ = 02/01 = 0.23 exp(— i32°) from Eq. 1)4. 6(1 . 
Evidently, the ratio of C/T in i3 — > vrvr decays is similar to that in B ^ Dtt decays [see Eq. (|3.43|) ]. 
To analyze the B tttt data, it is convenient to define %,g = T + £ + V and Cgg = C — £ — V . It 
follows from Eq. ()5.17|) that 



Ceff 



0.71 e 



j72° 



(5.18) 



This is consistent with Eq. (|5.9j) obtained from the global analysis oi B ^ tttt and vri^ data. 



C. CP violation in i? ^ vr vr° decays 

We see from Table |IV| that for p- TT iT^, even including FSIs CP- violating partial-rate 
asymmetry in the SM is likely to be very small; the SD contribution ~ 5 x 10~^ may increase to 
the level of one percent. Although this SM CP asymmetry is so small that it is difficult to measure 
experimentally, it provides a nice place for detecting New Physics. This is because the isospin of 
the vr~vr'' ( also p^) state in the B decay is / = 2 and hence it does not receive QCD penguin 
contributions and receives only the loop contributions from electroweak penguins. Since these 
decays are tree dominated, SM predicts an almost null CP asymmetry. Final-state rescattering 
can enhance the CP- violating effect at most to one percent level. Hence, a measurement of direct 
CP violation in B~ vr~vr'^ provides a nice test of the Standard Model and New Physics. This 
should be doable experimentally in the near future. If the measured partial rate asymmetry turns 
out to be larger than, say, 2%, this will mostly likely imply some New Physics beyond the Standard 
Model. 
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VI. B ^ pTT DECAYS 



The experimental results of CP averaged branching ratios and CP asymmetries for B ^ pn 
decays are summarized in Table [V] The experimental determination of direct CP asymmetries for 



p~^Tr and p vr^ is more complicated as B^ 
CP asymmetries are given by 



p^TT^ is not a CP eigenstate. The time-dependent 



Y(B^(t) 



{Sp-n ± A5p7r) sin(Amf) — (Cpyr ± ACpyr) cos(Amt), 



(6.1) 



where Am is the mass difference of the two neutral B eigenstates, Sp-p^ is referred to as mixing- 
induced CP asymmetry and Cp^r is the direct CP asymmetry, while AS'pvr and ACp^r are CP- 
conserving quantities. Next consider the time- and flavor- integrated charge asymmetry 

iV(p+7r-) - N{p-T:+) 



N{p+Tl-) + N{p-TT+) 

where N{p'^7r~) and N{p~7r^) are the sum of of yields for B^ and 5^ decays to p'^'K^ and p^vr 
respectively. Then, 

T{B^ 



(6.2) 
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TiB" 
r(5° 



p+TT-)+T{B^ p-7^^ 

p^7T+) - r{B° p+TT- 



A 



c. 



1- AC, 



*^p7r^p7r 
H~ *^p7r^Cp7r 



p TT"' 



T{B p-7r+) + r(50 ^ p+TT- 



1 -|- ACpji- -|- ^p-TfCpTT 



(6.3) 



Note that the quantities ^p±7r=F here correspond to defined in 71 1. Hence, direct CP asym- 



and „4p-7r+ are determined from the above equations together with the measured 



metrics .4p+^- 

correlation coefficients between the parameters in the time-dependent fit to B^ 
Table O it is evident that the combined BaBar and Belle measurements of B^ - 
3.6(7 direct CP violation in the p^vr" mode, but not yet in p~tt~^ . 
The B ^ pn decay amplitudes have the general expressions 
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p^vr^ imply a 



AiB^- 


> p^vr ) 


A(^- 


* p~7r^) 


a{b'- 


- pV°) 


A{B- - 


- p-vr") 


A{B~ - 


- P°vr-) 



1 



1 



{Cv-Vv + r^'^ + ^r'^'^ + Cp 



Vp + vf^ + 



1 



-j={rp + Cv + Vp 



Py + Py 3^ 3^ 



V2 



(Ty + Cp + 



Pp + Vf^ + 



V 



They satisfy the strong isospin relation 



V2 



A{B- 



p vr 



+ A{B- ^ p^t:~) + V2A{B^ 



AiB" ^ p+TT~) + AiB^- 



P TT 



(6.5) 
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TABLE V: Experimental data for CP averaged branching ratios (top, in units of 10^®) and direct 
CP asymmetries (bottom) for B ^ pir 
correspond to defined in [t^I- 



50|,|65|. Note that the direct CP-violating quantities A 



Mode 


BaBar 


Belle 


CLEO 


Average 




22.6 ± 1.8 ±2.2 


29.lt|^±4.0 


27.6l^:^ ± 4.2 


24.0 ±2.5 




1.4 ±0.6 ±0.3 < 2.9 


5.1 ± 1.6 ±0.8 


< 5.5 


1.9 ± 1.2'* 


B- p-Tr° 


10.9 ± 1.9 ± 1.9 


13.2±2.3tJ:^ 


< 43 


12.0 ±2.0 


B- pOyr- 


9.5 ±1.1 ±0.9 


8.0t||j ±0.7 


10.4t||±2.1 


9.1 ± 1.3 


Ap~-^+ 


-0.21 ±0.11 ±0.04 


-0.02±0.16l[]:[J^ 




-0.15 ±0.09 


Ap+T^- 


-0.47 ±0.15 ±0.06 


-0.53 ± 0.29t[};[j^ 




-0.47l°:J3 


Ap-j^o 


0.24 ±0.16 ±0.06 


o.o6±o.i9j:|]:[|^ 




0.16 ±0.13 


ApO^j- 


-0.19 ±0.11 ±0.02 






-0.19 ±0.11 



'^The error-bar of B{B p^n ) is scaled by an S factor of 1.9. We have taken into account the BaBar 
measurement to obtain the weighted av erag e of the branching ratio for this mode. Note that the average for 
p°7r") is quoted to be < 2.9 in 



In the factorization approach (the subscript "SD" being dropped for convenience) 



Typ = KvpXuO.!, 



Cyp = Ky p\ua2: 'Pp = l^P 



Vy = Ky 



-n'EW 
' V 



^"(as^x + flio) + ^c(a8^x + ^^o^ 



' V,P 



3 



K'V,p{K + Ac)(-a7 + ag) 



-Kp 
2 



10 



with 



ny = V2GFf7zAQ''{ml)mp{e ■ pp), Kp = V2GFfpFf''{ml)mp{e ■ pp), 
Short-distance contributions yield 



B(B^ - 


> P vr+)sD 


= 18.4 X lO^*', 


A^^ ^ 


= -0.03 


B(B^ - 


> /0^vr")sD 


= 7.9 X 10-^ 




= -0.01, 


B(B^ - 


- P°vr°)sD 


= 5.9 X 10~^ 




= 0.01, 


B{B' - 


/0"vr°)sD 


= 12.8 X 10~^ 




= -0.04, 


B{B- - 


/vr~)sD 


= 6.8 X 10"^ 




= 0.06 . 



(6.6) 
(6.7) 



(6.8) 

Note that the central values of the branching ratios for -n^ , p'^ t[~ , p^jr^ , p~tt^, p^tt~ , respectively, 
are predicted to be 21.2, 15.4, 0.4, 14.0, 11.9 in units of 10~^ in [l^, which are larger than our 
results by around a factor of 2 for p^vr" and 1.5 for p^n~ . This is ascribed to the fact that the form 



factors AQ^iO) and Ff^'iO) are 0.28 and 0.25, respectively, in our case |63], while they are chosen 



to be 0.37 ± 0.06 and 0.28 ± 0.05 in fl3|. The central value of the predicted B(^ 



(36.5 



+21A-, 
-17.6) 



X 10 [l^ is large compared to the experimental value of (24.0 ± 2.5) x 10 , but it is 



accompanied by large errors. 
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(a) 



(b) 



FIG. 8: Long-distance t-channel rescattering contributions to -B — > pvr. 



The long-distance contributions are depicted in Fig. |H| 
Gf . \pi\dcose 2 F^{t,m^) 



V2 J-i lovrms 



X -P2 • P3 + 



(pi •P2)(P1 -Ps)^ ^2|P4| - E4\p2\ COS 9 



'af^A^,P{ml) + hf,F^-{ml) 



uibIpa 



Abs{8b) = -j=Xu / —r^ 2\/2 qd-Dtt 9v P — — 

V2 J-i IQirniB t — mj^ 

X 'cfDA^''*{ml)+dfD*Ff''{ml.) 

{pi ■P2){P1 ■P3)\ E2\P4\ - E4\P2\C0S9 



X -P2 ■ P3 + 



mi 



Abs{8c) = ^ AV2gl,DngvX 



mBlPA] 
t - mjj. 



-mo* 



cfDA^^*{ml)+dfD*Ff^{ml,)\ [mlpi ■ p2 - {P2 ■ Pa){pi ■ Pa) 



+ 



El \pa I + -£'4 \pi I cos 6 



rriBlPAl 



where a = 6 = 1/2 for p^ir^ 



(PB ■ P2){P3 ■ Pa) - {PB ■ Pa){P2 ■ Pd] 



(6.9) 



1(1/^/2), 6 = for /3^7r+ (/3°7r-), 



0, b = l{l/V2) for 



p+vr" (p-7r°), c = d= 1/2 for p^Tr^, c = d= 1/V2 for p^-K' , c = d = -1/V2 for p"7r°, c = l,d = 
for p~iT^, and c = 0, d = 1 for p~^7r~. 

A fit of the cutoff scales A7r,A^(*) or the cutoff parameters 7]^^ and ri^(,)to the measured 
branching ratios of — > pvr shows a flat without a minimum for r]£, < 2. Therefore, we choose 
the criterion that FSI contributions should accommodate the data of the p^ir^ and p^ir^ without 
affecting p^vr^ and p^vr*^ significantly. Using rjD = rjD* = 1.6 , the branching ratios and direct 
CP asymmetries after taking into account long-distance rescattering effects are shown in Table IVTl 
(The results are rather insensitive to 77,^.) We see that, in contrast to the ttK case, the B ^ pn 
decays are less sensitive to the cutoff scales. The B rate is enhanced via rescattering 

by a factor of 2, which is consistent with the weighted average of (1.9 ± 1.2) X 10-6. However, it 
is important to clarify the discrepancy between BaBar and Belle measurements for this mode. It 
should be remarked that direct CP violation in this mode is significantly enhanced by FSI from 
around 1% to 60%. Naively, this asymmetry ought to become accessible at B factories with 10^ 
BB pairs. As stressed before, there is 3. 60" direct CP violation in the p~^7r~ mode. Our prediction 
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TABLE VI: Same as Table HIl except for B ^ pn decays. 



Mode 



Expt. 



SD 



SD+LD 



B(B" ^ p-7T+) 


13.9+^-? 


18.4 


i8.8i°:i 




10.11?:^ 


7.9 


8.4 ±0.3 


B(B° ^ p^TT^) 


24.0 ±2.5 


26.3 


27.3in 




1.9 ±1.2 


0.6 


1 q+0.4 

-"-•-^-o.s 


B{B- p-vr") 


12.0 ±2.0 


12.9 


i4.ol°:I 


B{B- p°7r~) 


9.1 ± 1.3 


6.8 


7 c:+0.6 
' -^-0.3 




-0.15 ±0.09 


-0.03 


-0.24 ±0.6 




n 47+0-13 
"•^'-0.14 


-0.01 


-0.43 ±0.11 






0.01 


U-'^' -0.03 


Ap-j^o 


0.16 ±0.13 


-0.04 


0.36 ±0.10 




-0.19 ±0.11 


0.06 





of Ap+^- agrees with the data in both magnitude and sign, while the QCD factorization prediction 
(0.6l^Jj^ g)% |l^ seems not consistent with experiment. 



VII. POLARIZATION IN B ^ (pK*, pK* DECAYS 



The branching ratios and polarization fractions of charmless B VV decays have been mea- 
sured for pp, pK* and (pK* final states. In general, it is expected that they are dominated by the 
longitudinal polarization states and respect the scaling law, namely 2^ (see also footnote 1 on p. 3) 



1 



fL = 0{ml/ml), fjf\\ = 1 ± 0{mv/mB) 



(7.1) 



with fL,f± and /y being the longitudinal, perpendicular, and parallel polarization fractions, re- 
spectively. However, in sharp contrast to the pp case, the large fraction of transverse polarization 
in B — > (j)K* decays observed by BaBar and confirmed by Belle (see Table rvTl|) is a surprise and 
poses an interesting challenge for any theoretical interpretation. The aforementioned scaling law 
remains true even when nonfactorizable graphs are included in QCD factorization Therefore, 
in order to obtain a large transverse polarization in i? ^ 4>K* , this scaling law valid at short- 
distance interactions must be circumvented in one way or another. One way is the Kagan's recent 
suggestion of sizable penguin-induced annihilation contributions 23]. Another possibility is that 
perhaps a sizable transverse polarization can be achieved via final-state rescattering. To be more 
specific, the following two FSI processes are potentially important. First, B — > D*D* followed by a 
rescattering to 4>K* is not subject to the scaling law and hence its large transverse polarization can 
be conveyed to (pK* via FSI. Second, the FSI effect from B D*Ds or i? — > DD^ will contribute 
only to the A± amplitude. In this section, we shall carefully examine the long-distance rescattering 
effects on (j)K* polarizations. 
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TABLE VII: Experimental data for CP averaged branching ratios (in units of 10~^) and polariza- 
tion fractions for B —>■ (f)K* , pK* , pp [sO, 65 , [tJ, • The phases of the parallel and perpendicular 
amplitudes relative to the longitudinal one for B (j)K* are also included 0) • 



Mode 


BaBar 


Belle 


CLEO 




Average 




9.2 ±0.9 ±0.5 


10 n+i-6+0-7 




1.8 
1.7 


9.5 ±0.9 


B- (j)K*- 


12.7l:li ± 1-1 


7+2.1+0.7 
0- '-1.9-1.0 




1.8 
1.6 


9.7 ± 1.5 


B° p+K*- 


l6.3±5A±2.3tli < 24 








< 24 


B- p^K*- 


10.6l|^ ± 2.4 




< 74 




10.610 


B- ^ 


17.0±2.9±2.0j:?:^ 


6.6 ±2.2 ±0.8 






9.6 ±4.7'^ 


B- p^p- 


22.5t|:J ± 5.8 


31.7±7.ll|? 






26.4111 


B -^p+p 


30 ± 4 ± 5 








30 ±6 


B^ p^uj 


12.6l|^ ± 1.8 








12.611^ 




0.52 ±0.05 ±0.02 


0.52 ±0.07 ±0.05 






0.52 ±0.05 




0.22 ±0.05 ±0.02 


0.30 ±0.07 ±0.03 






0.25 ±0.04 


fL{4>K*-) 


0.46 ±0.12 ±0.03 


0.49 ±0.13 ±0.05 






0.47 ± 0.09 


f±{4>K*~) 




0.12l°;J^ ±0.03 






12+°-^^ 

^•-"-^-0.09 


(p\\{(l)K*^) (lad)^ 


2.34t[}:^ij ± 0.05 


-2.30 ±0.28 ±0.04 










2.47 ±0.25 ±0.05 


0.64 ±0.26 ±0.05 








(t)\\{^K*+) (rad)'' 




-2.07 ±0.34 ±0.07 






-2.07 ±0.35 


(l)±{(l)K*+) (rad)'' 




0.93t|]:i ±0.12 










0.96l|]:?^ ± 0.04 












0.79 ±0.08 ±0.04 ±0.02 


0.50±0.19i[]:|]f 






0.74 ± 0.08 


fLip'p-) 


o.97l:[5:[j? ± 0.04 


0.95 ±0.11 ±0.02 






0.96in^ 


hip^p-) 


0.99 ±0.03 ±0.04 








0.99 ± 0.05 


fUp'^) 


0.88lE]:J^ ± 0.03 











-*0, 



"The error in B{B- -> p-K ) is scaled by S = 2.4. Our average is slightly different from the value of 
9.2 ±2.0 quoted in 0. 
^Experimental results of the relative strong phases are for (j}K*'^^'^'> decays. 



A. Short-distance contributions to i? ^ (pK 



In the factorization approach, the amplitude of the B — > (f)K* decay is given by 

g= 

\BK* 1^2- 



{4>{X^)K*{\k*)\H,^\B) = E mss)v\Q){K*\{sh)v-A\B) 

n=u,c 



[m-j 



{niB + niK* 



jPBfiPBu 



2V'''' (ml) 



{vfiB + rriK' 



a 13 
y2^lJ,vapPBPK 



(7.2) 
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where A^^i^-* are the corresponding hehcities, a^^^, = a^i + a\ + — (07 + 09 + afo)/2 and 
a^K* = Gprn^f^imB + mK')y.„-,. . VnhV*. a^^,/\/2. With Ax defined as {4>{X)K* {\)\Hcs\B) , 
it is straightforward to obtain [2a] 

A± = -ia^K' Af"^ {ml) ^ ^' niBPc , (7.3) 

where pc is the center of mass momentum and we have used the phase convention e^{X) = £k*{—X) 
for p^ji^K* — *■ 0. The longitudinal amplitude ^0 is sometimes denoted as A^, while the transverse 
amplitudes are defined by 

A\\ = = -ia^K'^A^^' 

in the transversity basis. The decay rate can be expressed in terms of these amplitudes as 

- -^{\Ao\' + \A+\' + \A.\') = -^{\Al\' + l^iil^ + \A^n (7.5) 



and the polarization fractions are defined as 
with a = L, II, _L. 

The form factors A,.A,.V usually are similar irr size 0. Using above equatiorrs with a. 
given in 1)4. 6|) and the form factors given in [g^ ]. one immediately finds that the short-distance 
contributions yield (see also Table rvTTl|) 

/!° : /|f° : /1° = 0.88 ± 0.06 : 0.07 ± 0.03 : 0.05 ± 0.03, 1 ~ /|° > (7.7) 

where the errors are estimated with 10% uncertainties in form factors. However, the above pre- 
diction is not borne out by the B — > (pK* data (see Table IVIip . It is easily seen that Eq. ()7.7() is 
closely related to the fact that ms ^ m^^K*- We do not expect that such a relation holds in the 
case of — > D^*^D^s^ decays, which can rescatter to the (pK* final state. It is thus interesting to 
investigate the effect of FSI on these polarization fractions. A similar analysis was also considered 
recently in j?^ . 

Another possible shortcoming of the short-distance factorization approach is that the central 
value of the predicted branching ratio 

B(B ^K*)sD = (4.4 ± 2.6) x 10~^ (7.8) 



As shown explicitly in [t^ within the QCD factorization approach, nonfactorizable corrections to each 
partial-wave or helicity amplitude are in general not the same; the effective parameters a.i can vary for 
different heli city amplitudes. However, since such an effect is too small to account for the large transverse 
polarization (26j , for simplicity we shall ignore the difference of Oj in different helicity amplitudes. 
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FIG. 9: Long-distance t-channel rescattering contributions to i? ^ (pK* . Note that flows in lines 
are along b or c quarks. 



is too small by a factor of 2 compared to experiment. As pointed out in 77], in the heavy quark 
limit, both vector mesons in the charmless B — > VV decays should have zero helicity and the 
corresponding amplitude is governed by the form factor ^o- Consequently, the predicted branching 
ratios are rather sensitive to the form factor models of Aq and can be easily different by a factor 
of 2. Using the covariant light-front model for B K* form factors |6^, we are led to the above 
prediction in Eq. 1)7. 8|) . A similar result is also obtained in the BSW model [t^I. However, a much 



7£,B) 



larger branching ratio will be obtained if the form factors based on light-cone sum rules 
are employed. At any rate, if the predicted rate is too small, this is another important incentive 
for considering long-distance effects. 



B. Long-distance contributions to B ^ (j)K 

The states from B decays can rescatter to (pK* through the t-channel -Oi*^ exchange in 

the triangle diagrams depicted in Fig. IHl The effective Lagrangian for £)(*)£)(*) y vertices IS given 
in Eq. 1)3.15(1 . Note that as in Eq. (|4.12|) . SU(3) or U(3) breaking effects in strong couplings are 
assumed to be taken into account by the relations 



9d(*)D^*^K* = \ -Z~SD(*)Di*)p: = " QDMDMp^ (7-9) 

and likewise for and /^(.)^(.)^. 

There are totally eight different FSI diagrams in Fig. The absorptive part contribution of 
B DgD — > (j)K* amplitude via the Dg exchange is given by 

MsiD.D-M = \ j ^^^^11^^ ^^^11^^ i2.mps - - pM(B^ D.D) 

xm9DsDs<^^Y^'^{-2i)9DDsK* {el-Vi){el-P2). (7.10) 

We need to use some identities to recast the above amplitude into the standard expression given in 
Eq. (|7.2() . To proceed, we note that after performing the integration, the integrals 

/ (2^)5^1 (2^)S^2 (^^^^'^^'^^^ -Pi- P2)f{t) X {pi^, Pli^Plu, PlpPluPla} (7.11) 
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can be expressed only in terms of the external momenta P3,Pa with suitable Lorentz and permutation 
structures. Therefore, under the integration we have 

PltjiPluPla = {QfjivPa + QfiaPy + gvaPfi)^^^ + {gf^uQa + Qfiaqv + guaqfi)^^^ 

+P^P,Po,Af + {P^P.q,, + P^q.Pa + q^PuPc.)Af 

+ {qfj.quPa + q^iPuqa + Pf^qvqa)Af^ + q^qyqaAf\ (7.12) 

where P = p^ + P4, q = P3 — Pi and A^j'^ = A^^ {t,m\,m\,m'^,m'l,rn\) have the analytic forms 
given in Appendix A. With the aid of Eq. 1)7. 12(1 . we finally obtain 

AbsiD.D;D.) = \ J ^^^tl^^ ^^f^^^ i2nmPB - P. - pM(^ ^ D.D) 

F\t,mDj 

X4 gDsDs<f>—. —2 9DDsK* 

^ ''''Ds 

x{el . eli-Af^) - {el • P){el ■ P)(-^S^) + + - )}. (7.13) 
Likewise, the absorptive part of the B DgD — > (j)K* amplitude via the D* exchange is given by 

15 

Abs {DsD; D*) = - f , "^y^ , (2^)^<5^(PB - pi - P2)A{^ ^ D,D) 

F2(t,mcj) 

X lb/D,DJ0— —2 JDD*K* 

i III r)* 



-e^ -Pel-P 



2A?^-p,.pM'^'-A?-A?^+Af^)]\, (7.14) 



where k = pi — p^. It can be easily seen that the FSI contribution from the B DDg decay will 
affect both Al and A\\ amplitudes of the B (pK* decay. 

The absorptive part contributions of i? ^ D*D (f)K* amplitudes via Dg and D* exchanges 
are given by 

Ah m*n ^ [ ^Pi^P2,^ .4.4. A^" ^ DID) 

Abs[DgD;Ds) = - / ,^ .^^^ ,^ (27r) 6 [pb - Pi - P2) 



2 J (27r)32Ei (27r)32^2 ^ " 2e* • P 



To avoid using too many dummy indices, we define [A, B,C, D] = ea0jsA°'B^C'' , [A, B,C, ^] 
f-ap-ifiA" B^C^ and so on for our convenience. 
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(27r) 6 [pB -pi -P2)- 



Abs{DlD;Dt) = 



2 J {2n)^2Ei {2'!Tf2E2 



2el P 



t — mj) 



2 9DDsK* A[ 



2 J {2-K f2Ei (27r)32^2 



(27r) 5 (jpB -P1-P2) 



A{B DID) 



2e\ P 



m 



1 /■ dVi c?V2 



/ — 
y (27r 



2i (27r)32£;i (27r)32S2 



(27r)V(pB-pi-j92) 



^(5^ ^ D*D) 
2£^ -P 



x(-i)[e^,£^,P3,P4] (-32i)/n.n* 



t — m|). 



+ (l-4'^-4'^)(P-P3 



I (2) 



c7 + (l-a 



P-pi 



2 



P • pipi • P3 



(7.15) 



where the dependence of the polarization vector in ^(P — D*D) has been extracted out exphcitly 

and cr = gD*D*v/i2fD*D*v)- 

Similarly, the absorptive parts of P — > DgD* —>■ cpK* amplitudes via Dg and D* exchanges are 
given by 



Abs {DsD*-Ds) 



If d^pi d^p2 .0^4.4. AiB^^D.D* 

(27r) 6 [PB-Pi -P2)- 



2 J (27r)32Pi (27r)32P2 



2e*.P 



^2^2 ■P2igD,D,4> I J ^2 " {-^i)fD*D,K* \P4:,£hP2,£2]pi ■ 



A2 



1 f d^pi d^p2 



A{B^ -> DsD* 



(27r)32Pi (27r)32P2 ' ' ' ' 2£^ P 

X (-i) [£* , £1, P3, P4] { 16i W 4'^ I , 



Abs {D,D*;D*) 



1 /■ d3pi d^P2 



I — 

J {2tt 



{2T:f5\pB-pi-p2) 



2 J (27r)32Pi (27r)32P2 



A(^ DsD*) 
2i*TP 



X ^ 2£^ P 4i fo^Di^ ^_^2 
A2 



-{'ii)fD*D*K* 



D* 



( k^k^\ 



1 /• (i3pi d3p2 .4e4. A(P° ^ P>,P>* 

(27r) 5 (pB -pi -P2)- 



27 (27r)32Pi (27r)32P2 



2£* • P 
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-fD'DtK^{A\ 



,(2) 



m. 



P ■ P2P2 ■ Pi 



ml 



2 

(7.16) 



where again the dependence of the polarization vector in A{B DgD*) has been extracted out 
explicitly. Note that through rescattering FSI, both B D*D and B — > DgD* will affect only the 
A_\_ term of the B — > (l)K* decay amplitude. 

To consider the FSI effect from the decay B ^ D*D*, we denote the decay amplitude as 

A(B ^ D:{pi,Xi)D*{p2,X2)) = £T4"{a9^.. + bP,,P, + ic[^, „ P,p2]). (7.17) 

The absorptive part of the B D^D* (pK* amplitude via a t-channel Dg exchange has the 
expression 

2 J (27r)32^i (27r)32£;2 ^"^^^^^^^P^ -P^- P2){ag„v + bP^,Pu + ic[^, P,P2]) 

F^{t,mDj t A-\^ \ * IF * 1 

{-^l)jD*DsK* m,£2,^Pl,£l\\PA,£i:P2,£2\ 



Ahs {DlD*-Ds 



J2 eTer i-moiD. 



Ai,A2 



t — m 



D, 



1 



(fpi d?p2 , .44, . F'^{t,mDj 

(27r) 6 [PB - Pi - P2)lQfD*D,(t>— — fD'D, 



2 J (27r)32£;i (27r)32£;2 



K* 



X |e* • el{2aps ■ p.A^^- [{ps ■ p^f- mimi] [a{A\'>- A\''- A),'>+ <0 + hA'^']} 

-el ■ Pel ■ P{2aAf'^ -p^ ■ Pi[a{A^^^ - - A^? + 4'^) + ^'^i'^]} 

-i[el, el,p^,p^] 2c{P • P3(4'^ - 4'^ -A§^)+P- pM^^ - 4'^) - Pz ■ PiA^2^ 



2_2ir / .(I) .(1) 4(2)^ /l(2)^ 



1(2)1 



+ 



[(P3 • Pi? - mlml]{Af^ - Af - Af + 4^^)} 



(7.18) 



We see that the above FSI contributes to all three polarization components. The contribution from 
the D* exchange is 



Abs {D*D*;D*) 



1 f d?pi (f'p2 



2 J (27r)32Ei (27r)32£;2 



(27r) 5 {pB-pi-P2){agsp + bPsPp + ic[s,f3,P,p2]) 



^i^^^'^ (.-^i)fD*Dt<i> ^. ^ {^i)fD*D%K* e{{agp^pi ■ el+p^pel^- elpP^^) 



Ai,A2 



X ( - 5'^" + 2 



{(^9uaP2 ■ el - PivSlo, + £3^P3a)£? 

1 [ Spi d?p2 ,^ ,4^4, F'^{t,mDx) 



1 [ (Ppi (Pp2 , .4„4. . , 

y I L rrijj* 



, _ 2 fD*D*K* 



aP3 -Pi + bps ■P2P1 -PA 



a + bpi-p2, 2 2 , ^' 

H 0-^5 — (-"^2^3 • PiPi ■ Pa - m^ps ■ P2P2 ■ P4 + Pa • PiPi ■ P2P2 ■ Pa) 



2 2 

m J 771,2 



+a -P3-PA + 



P3- kk- Pi 
ml. 



+ AfU 



o u a + bpi-p2f .2 
2a - bpi ■ p2 H 2~^ — [pi ■ P2) 

mfm2 
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a + hpi ■ p2 
b 2 2 Pi • ^'2 

mfm2 



. a + ^'Pi • P2 . 2 2 

0P2 ■ P3 H 2~~2 ("^iK • P3 - m2P\ ■P3-P2- PlPl " Ps) 

777/ 777/2 

^ a + • P2 . 2 2 , 

- bpi ■ P4 H 2 2 ("^lP2 • P4 - ^2^1 -PA+Pl- P2P2 ' Pa) 



+ P3 • P4 - 



Ps- kk ■ p4 



m, 



Dt 



b + 



mfmg 

a + bpi ^P2 ^^2 _^^2_^p^ .^2) 



l"''2 



a + 6pi • p2 



(1-a) 



m 



"(4 



(1) 



(1) 



ap3 ■PA + bp3-p2Pi - Pa 



2 2 



(-m|p3 • pipi • p4 - mfp3 • P2P2 -PA + Ps- PiPi ■ P2P2 ■ Pa) 



(2) 



2a - bpi ■ p2 + ^^^__2^'^l {a + bpi ■ P2) 



m\m2 



■ o + bpi ■ P2 

b — pi ■ p2 

mfm2 



-(4'^-AW)(P1-P4+P2-P3)]- 

a + bpi- p2 



Dt 



a + 



'-{-miP2 ■P3+P2- PiPi ■ Pa) 



l + {l-a){l-AP-A),'>) 



.(l)^P3^ 



m 



a + bpi-p2 
a + bpi-p4 5 p2 • Pa 



rrin 



4^)+4^)-(i-a)(4^)-Ar)^ 



?'P2 • P3 + 



a + • P2 . 2 2 

— {miP2 ■ P3 - m2Pi ■P3-P2- PlPl ■ P3) 



m\m2 



+ 



4^)-4') + (l-a)^(4^)-A 



m 



(1) 
2 



4^4^^; 



+ ^3-^4- 



a + 6pi • p2 

2 2 
m|r?T,2 

P3- kk- Pi 



bpi ■ Pa - —--T-^ZT-^i'^iP'i • ~ ^Ipi -PA+Pi- P2P2 ■ Pa) 



m 



Dt 



{6(4'^ +4^) -4^) +4^)) 



+^±5^bi • P2(A« - 4^) - 4^^ + 4^^) - miii - 2Af^ + 4^) 



2 2 

mfm2 
1(2) 



(2)> 



m^(4^^ - 4'^)] \ \ - i[elel,P3,PA] 2c\ [p^-p^ 



P3 ■ kk-p4 \^^[i) 



m 



Dt 



2ct - (1 - a) 



q ■ k 



A\ 



(2) 



(7.19) 



with P = P3 + Pa, Q = P3 — Pa and k = pi — ps = p4 — p2 as before. Clearly the above rescattering 
amplitude contributes to all three polarization components >1l || ±- 
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C. Numerical results 



In order to perform a numerical study of the analytic results obtained in the previous subsection, 
we need to specify the short-distance A{B — > amplitudes. In the factorization approach, 

we have 

i^ycbVcsaifDsirriB - m^)FQ (m^J, 
^VcfeVc* ai/D*mD*Ff ^(m|,*)(2ej5, • ps), 
^VcbV*saifD,rnD'AQ^'' {m%^){2e*o* -pb), 
-i^VcbV*gaifD*rnD*{mB + mD*)e*j^.e% 



A{B - 


^D,D)sD = 


A{B- 


^DlD)sD = 


A(B - 


^DsD*)sD = 


A(B- 


>d:d*)sd = 



[mB+rnD*) {mB+mD*y 

Writing A = A^^ + iAbsA^^ with AbsA^^ being given in the previous subsection and using the 
form factors given in ,63;], we obtain the numerical results for B — > (j)K* amplitudes as exhibited 
in Table IVIIII The branching ratio after the inclusion of FSIs is 



2V^^ (ml,) 



a P 
7^tiuaf3PBPD- 



(7.20) 



B{B ct>K*) = (9.2l|^) X 10 



1-6 



(7.21) 



for A£)^^£)| = rriBi^^^D* +?? Aqcd with 7] = 0.80 and 15% error in Aqcd- Several remarks are in order, 
(i) The polarization amplitudes are very sensitive to the cutoffs Ajj^^d*. (ii) We have assumed a 
monopole behavior [i.e. n = 1 in Eq. (|3.22|) ] for the form factor F{t,m£)J and a dipole form 
(i.e. n = 2) for F(t,m£)»)}^ If the monopole form for F{t,mB)t) is utilized, the parameter rj 
will become unnaturally small in order to avoid too large FSI contributions, (iii) the interference 
between short-distance and long-distance contributions is mild. This is because B (pK* and 
B —>■ amplitudes should have the similar weak phase as |V"cfeV^*| ^ iKib^usI- ^ 

result, the long-distance contribution is essentially "orthogonal" to the short-distance one. (iv) 
The calculated strong phases shown in Table Ivnil are in agreement with experiment for cpii but not 



For B Dtt, Ktt, tttt decays, assuming a dipole form for the form factor appearing in the VcxcVP or 
VcxcVV vertex with Vcxc being the exchanged vector particle [see e.g. Figs. EIe),0If),El[d) andEl^e)] and a 
monopole behavior for all other form factors, we found that the best fit to the measured decay rates is 
very similar to the one with the monopole form for all the form factors. Only the cutoff A or the parameter 
r] appearing in Eq. H3.40|l is slightly modified. The resultant branching ratios and CP asymmetries remain 
almost the same. The same exercise also indicates that if the dipole behavior is assumed for all the form 
factors, then the value becomes too poor to be acceptable. This is the another good reason why we 
usually choose n = 1 in Eq. (|3.22() for the form factors except for the one appearing in VcxcVP or VcxcVV 
vertex. 
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The weak phase of SD B (j>K and B Dl* D^*^ is not exactly the same. Moreover, our B ipK 
SD amplitude is not purely imaginary. It carries a phase as shown in Table IVnil mainly because of the 
phase difference in and in our case [cf. Eq. (|4.6|l ] . 
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TABLE VIII: Short-distance (SD) and long-distance {LD) contributions to the decay amphtudes 
(in unit of 10^*^ GeV) and the polarization fractions for B — > (pK* . LD contributions are calculated 
with Ad^^d* = rriD^^D* + r/AqcD) where i] = 0.8 and errors come from a 15% variation on Aqcd- 





longitudinal (L) 


parallel ( ) 


perpendicular (±) 




0.63 - 1.98i 


-0.18±0.55i 


0.15 - 0.46i 


j^SD+LD 


-0.57toil - (2.03±0.01)i 


2.28t[!:^^ + (0.64 ± 0.03)i 


0.32 ±0.05 - (0.44 ±0.00) 


A{DsD;Ds) 


-l.lOtlil - (0.04±0.01)i 


0.07 ± 0.02 ± O.OOi 





A{D,D-Dl) 


0.00 


-0.02 ±0.01 - O.OOi 





A{D*D;D,) 








-0.19 ± 0.05 -O.Oli 


A{DID-Dl) 








-0.01 - O.OOi 


A{D,D*;Ds) 








0.17 ± 0.05 ± O.Oli 


A{DsD*;D*) 








0.01 ± O.OOi 


A{DtD*;Ds) 


-0.06 ±0.02 -O.OOi 


2.35l°:el + (0.09l°:°i)z 


0.211°:°^ ± O.Oli 


A{D*D*;Dl) 


-OMto^l - O.OOi 


0.06lo:o3 ± O.OOi 


-0.01 ±0.01 - O.OOi 


fSD 


0.88 


0.07 


0.05 


jSD+LD 




0.541°:!° 


0.03 ±0.01 


yexpt 


0.51 ±0.04 " 


0.27 ±0.06 " 


0.22 ±0.04 


arg(^MI)^^ 




vr 





aTg{AiAlf^+^^ 




9 19+0.11 

^•-■-^-o.o? 


u.oy_Q 24 


arg(^i^^)^'^'^'^'' 




2.34l°;^3 ± 0.05 


-0.67 ±0.25 ±0.05 ^ 


arg(A,^2)'''"' 




2.21 ±0.22 ±0.05 


-0.72 ±0.21 ±0.06 ^ 



Qjoxpt jg ^j^g experimental average of polarization fractions for ^ (f>K*^ and B (pK 
''See footnote 18 for details. 



for i;^_L.^^ (v) There are large cancellations among various long-distance contributions, e.g. between 
A{DID-D^:^) and A{DsD*-D^:^). 

The cancellation in ^ D^D — > (f)K* and B DgD* — > (pK* contributions can be understood 
from SU(3) symmetry, CP conjugation and the similarity in the size of source amplitudes. To 
see this, let {4>{ss)K* (qs)\Tst\D*{cs)D{qc)) be the D*D (f)K* rescattering matrix element, which 
is represented in the right-hand part of Fig. |^ After changing s ^ q m. D*{cs) and 0(ss) (on 



In Table IVni the experimental results of 0|| and 0x are for B (j)K* decays. The Belle results for ((/)|| , (j>±) 
should be transformed to i—<t>\\ , tt — (t>±)- This phase transformation will not modify the experimentally 
measured distributions (see [8l| for details) but it will flip the relative sign between (/)||, (j)± and lead to 
\A^\ > \A^\ for B — > (pK* as expected from the factorization approach. Since our calculations are for 
B (fK* decays, in Table IVlIll we have transformed BaBar and Belle results from (j)±_ to 0^ — vr so that 
\A+\ < \A^\ in B ^ <pK* . Note that in the absence of FSIs, 0|| = = n for B 4>K* and = tt, 
(f)±_ = ioT B ^ (fK* . 

The combined symmetryis reminiscent of CPS symmetry which is used extensively in kaon decay matrix 
elements on the lattice l82Ll83l. 
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the upper right part of Fig. and q ^ s in D(qc) and K*(qs) (on the lower right part of of 
the same figure) under the SU(3) symmetry transformation, the above matrix element becomes 
{K*{sq)(j){ss)\Tst\D*{cq)Ds{sc)). Applying charge conjugation after the SU(3) transformation, we 
are led to 

{cP{ss)K*{qs)\T,t\D:{cs)D{qc)) = -{K*{qs)cl){ss)\T,t\D*{qc)Ds{cs)), (7.22) 

where the overall negative sign arises from the transformation of vector particles under charge con- 
jugation. In the B D*D (pK* rescattering amplitude, the above matrix element is sandwiched 
by e*{pi, A) -ps from B —>■ D*D amplitudes and efj_upa£^{p3, X')e'^{p4, X")p^p1 to project out the A± 
term. These factors bring the initial and final states into p-wave configurations. Applying parity 
transformation will interchange the particle momentum in the B rest frame. Such an interchange 
does not yield a further sign flip as the initial and final states are both in the p-wave configuration. 
To be more specific, we have 

Pj2\D*{p,,X)D{p2)) e{piArVBp = Y.\D*{P2.-X)D[pi)) [-6^(p2, -A)p^], (7.23) 

A A 

where the identity £fj_{p,X) = —£'-''{— p,—X) with p^ = {E,p) = p2^ in the B rest frame for the 
mi = m2 case and pi <-> p2 under the parity transformation (P) on initial state particles have been 
used. For the parity transformation of the final state, we have 

P J2 \K*{p3,X')HPi,X")) e^.p^e^{p3,X')e''{pA,X")p'3Pl 

A',A" 

= I^*b4,-A')</'(P3,-A")) (^fiupa{-f£f^iP4,-X')e^ip3,-X")p4pP3a 

A', A" 

= -Y^ \K*ip^, X")cP{p3, A')) e^,,,E>^{p3, A')e^(P4, A")p^pI, (7.24) 

A', A" 

where similar identities and indices relabelling have been used. We thus have 

E W-^^fe,A')e'^(P4,A")pgpJ {<Pip3,X')K*{p4,X")\Tst\D:{pi,X)D{p2)) e*{pi,X) ■ pb 

A.A'A" 

= - E ^^^upae''{p3,X')e•'{p4,X")p'3pl {K*ip^,X")<l>{p3,X')\T,t\D*{p2,X)Dsipi)) e*{p2,X) ■ pb 

A,A'A" 

(7.25) 

from SU(3) [U(3)] symmetry and CP conjugation. The right-hand side of the above equation is 
proportional to the B DgD* — > (pK* rescattering matrix element with an additional negative 
sign, which is the source of cancellation. This cancellation turns out to be quite effective as 
the B D*D and B DgD* amplitudes are numerically similar, i.e. fo'tno* 
fDsiTT'D*AQ^*{m'jjJ. Note that SU(3) breaking does not really help avoid the cancellation. Since 
the number of s quarks involved in B ^ D*D —>■ (f)K* and B — > DgD* (j)K* is the same, SU(3) 
breaking affects both amplitudes in a similar manner and the net effect cancels. 

Because of the above-mentioned cancellation and the fact that B D*D* does not give rise to 
a large A±, we cannot obtain large f_\_{B — > 4>K*) through FSI. The latter point can be seen from 
the fact that the polarization fraction of S ^ D*D* is /l : /n : /j_ ~ 0.51 : 0.41 : 0.08. Although 
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the D*D* intermediate state has a smaller Ai, it does not have a large A±. The suppression in 
A± can be understood from the smallness of the Pc/^B factor as shown in Eq. ()7.4() . In principle, 
one could obtain a larger Aj^{B i'K*) by applying a larger A/), to get a larger A{DID*] D*). 
However, such a large A{D*D*; D*) will also enhance Aj^^, giving a too large B — > (f>K* decay rate 
and does not really bring up f± . 



D. B- K*-p^ decays 

In the factorization approach, the amplitude of the B K*^ decay is given by 



{K*-{XK*)p''iXp)\H,fi\B- 



-| E ^ibV:, al.^ {p\{uu)y\Q){K*-\{sh)y.A\B~) 

q=u,c 

+% E y^^^K {K*-\{-su)v\0){p'\{ub)v^A\B' 

Q=u,c 



ax 



A^^\ml) 



{K*\XK*)p-{Xp)\H,^\B- 



-2i 



BK* (^2 ^ 
aK*p-, ; + Pk*p 



PBp.PBu 



+ (3k* 

{mB + mK*Y ''{niB+mp) 



(^puapP%Pp 



E ^'^bV:, C\.^ {r\sd)y\^){p~my.A\B- 
q=u,c 



{tub + rup)'^ 



PBpPBu - 2i 



m 



K* 



{irLB + nip 



a P 
^fiuafSPBPp 



(7.26) 



where 



'■K*p 



3 3 

02 + -(07 + 09), «K*p = 2(^7 + ag), bK*p = ai + a2 + a'lQ, 65^.^ = 0^ + 0^0, 



1 



1 



aK*p = ^fprnpimB+mK") E ^gbKjs "iC-p' Pk*p = ^fK^mK-'imB + nip) E ^qbVgs b^^f 



q=u,c 



q=u,c 



Gp 



7K*p = -j=fK'niK'{mB+nip) E ^9b^g*s c|-.p. 
Likewise, we have 



Ao{K*-p^^ 



[aK^pA^""* {ml) + OK'pA^^nil,)]- 



— mj^t — nip 
2mK*nip 



+2 



aK-p- 



Ar\ 



m 



A^Piml 



h PK*p-, ■ 72 

(ms + nipy 



Au{K*~p^) 



[niB + mx*)^ 
iV2[aK*pAf''\nil)+PK*pA^P{mh)] 



1 1 
rriBPc 

niK^nip 
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TABLE IX: Same as Table IVTllI except for B —>■ K* and Ad^jj* = "t-d.d* + ??Aqcd- 





longitudinal (L) 


parallel ( ) 


perpendicular (±) 




1.18 - 1.69i 


-0.26 ± 0.37i 


0.16 - O.lOi 


^SD+LD 


0.1ll°:i - (1.73l°;°?)z 


1.43t[}:^^ + (0-43 ± 0.02)i 


0.25 ± 0.03 -O.lOi 


A{D,D-D) 


-l.Ol^jjJJ - (0.04±0.01)i 


0.04 ± 0.01 + O.OOi 





A{DsD;D*) 


0.00 


-0.01 ±0.01 - O.OOi 





A{DID-D) 








-0.13t|]:[]^ - O.OOi 


A{D*D;D*) 








0.00 


A{D,D*-D) 








0.11 ± 0.03 ± O.OOi 


A{DsD*;D*) 








0.00 


A{D*D*;D) 


-0.03 ±0.01 - O.OOi 


1.62+°;!^ ± (0.06 ± 0.02)i 


0.12+°:°^ ± O.OOi 


A{D*D*-D*) 


-0.04 ±0.02 - O.OOi 


0.04lg;g^ ± O.OOi 


-o.oit[!:[j? - O.OOi 


fSD 


0.95 


0.04 


0.01 


j!SD+LD 
jex.pt 


0.57l°:J^ 
0.96lH^ 




0.01 



A^{K*~p^) = 2^/2i 



OiK*p-, ; ± Pk''p- 



rriBPc 



{mB + mK*Y {mB + rnpy_ 
Aq{K p ) = -i-/K*p[ -A/(m^*) ^±- , 

A\\{K*%-) = -t^K*pV2A^'{ml,), A^(K*%') = i^K'p '^^^Z^'^\^ ^^P^' (7-27) 
Similar to the calculation for B (pK* , we obtain 

B^'^iB- K*-p°) = (4.1 ± 2.2) X 10-^ 
e^^(S^^K*V) = (4.3 ± 2.9) X 10-^ 
/|° : /|f° : = 0.95 ± 0.02 : 0.04 ± 0.02 : 0.01 ± 0.00, 

/!° : /|f° : fl^(K*^P~) = 0.92 ± 0.05 : 0.04 ± 0.02 : 0.04 ± 0.01, (7.28) 
where the errors are estimated from 10% uncertainties in form factors. 

The formulae for long-distance contributions to K*~ p^ decays are similar to that in 

B (j)K* with the intermediate state D^*^ and strong couplings f {g) ^{*) ^{,) ^ and f{g)^^^^^(,)^^ 

replaced by D^*\ /(5r)^(.)^(,j^^ and f{g)oMD(*)p/^^ respectively. Likewise, the B' K*^ p' 
rescattering amplitudes can be obtained by a similar replacement except for an additional factor 
of l/^/2 associated with f (g) £)(*) d(') p- 

Numerical results for B^ ^ K*^ p^ and B^ ^ K*^p^ are shown in Tables HXl and 1x1 respec- 
tively. Using A^(,) = ± 7?Aqcd with -q = 0.8, we obtain 

^5D+LD(^*-^0) ^ (4_g+1.8) ^ 10-6^ e^^+^^(Z*V-) = (10.311?) X IQ-^, 

A^^{K*'p'^) = 0.30, A^'^iK*^^) = 0.11, 

^5D+LD(^*-^0) ^ o.02t°:|]^, ^^^+^^(K* V) - 0.00 . (7.29) 
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TABLE X: Same as Table IVmi except for K p~ . 





longitudinal (L) 


parallel ( ) 


perpendicular (±) 


^SD 


0.43 - 2Mi 


-0.10 + 0.45i 


0.09 - 0.41i 


j^SD+LD 


-imtlil - (2.09±0.02)i 


2.29t°l + (0.54tHi)^ 


0.22 ± 0.04 -0.41i 


A{D,D-D) 


-lA2to[^ - (0.05 ±0.02)i 


0.06 ± 0.02 + 0.00« 





A{D,D;D*) 


0.00 


-0.02 ±0.01 - O.OOi 





A{D*D;D) 








-0.18 ± 0.05 - O.OOi 


A{D*D;D*) 








-0.01 - O.OOi 


A{D,D*-D) 








0.16+g:|]| ± O.OOi 


A{DsD*;D*) 








0.01 ± O.OOi 


A{D*D*;D) 


-0.04 ± 0.01 -O.OOi 


2.29t°:™±(0.09l°:°^)^ 


0.171°:°^ ±o.oii 


A{D*D*;D*) 


-0.05l[J;[!| - O.OOi 


0-06to:o3 ± O.OOi 


-0.01 ± 0.01 - O.OOi 


fSD 


0.92 


0.04 


0.04 


jSD+LD 


4q+0-ii 

*-'-*^-0.08 


n 40+0.08 


0.02 ± 0.01 


jcxpt 


0.74 ±0.08 







Note that the long-distance contributions are similar to that of the (f)K* case up to some SU(3) 
breaking and the LD contributions to K*~ and K*^ p~ are related through the above-mentioned 
factor of from the D^*^D^*'^ p vertices. Before the inclusion of FSIs, K*^ p^ and 'K*^ have 
similar rates as the former receives additional tree contributions from T and C topologies. However, 
the LD contribution to the former is suppressed at the amplitude level by a factor of ^/2 and this 
accounts for the disparity in rates between these two modes. The cancellation between B PV 
VV and B — > VP VV amplitudes remains. Since the short distance B~ K*~ p^ amplitude 
has a non-vanishing weak phase, we have interferences between SD and LD contributions resulting 
in quite sensitive direct CP. Comparing to the data shown in Table IVIll we note that (i) the K*~ p^ 
and K p~ rates are consistent with the data within errors and (ii) the predicted fi {K p~ ) agrees 
well with experiment, while fL{K*^p^) is lower than the BaBar data. This should be clarified by 
future improved experiments. 

E. FSI contributions from other channels 

For LD contributions to the B — > (pK^*^ decay, it is also possible to have B J/ipK^*^ — > (pK* 
via a t-channel r/s (a pseudoscalar with the ss quark content) exchange and B K^*^J/%Ij cj)K* 
via a t-channel K exchange. However, it does not seem that it can fully explain the polarization 
anomaly for two reasons: (i) The CP and SU(3) [CPS] symmetry argument also leads to the 
cancellation in the B — > J/ipK — > (pK* via a t-channel r/^ exchange and B — > KJ/ip — > (j)K* 
via a t-channel K exchange contribution, (ii) The B — > J/ipK* decay does not generate a large 
enough A±_ source. The polarization fractions of ^ J/'ifjK* are fi '■ f\\ ■ f± — 0.47 : 0.35 : 0.17, 
which differ not much from that of -B ^ D*D*. On the other hand, it is known that SU(3) 
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breaking and the electromagnetic interaction take place in J /tp decays (in particular, in the case of 
J/ifj ^ K* K, K*K and (fyq^'^ decays of interest) at the level of ~ 20% and ~ 10% jsj], respectively. 
The above-mentioned cancellation may not be very effective. Hence, we should look at the FSI 
contributions from this channel more carefully. 

The Lagrangian relevant for J/ij: ^ PV and (p — > PP decays is given by 

C = igvPpTiiV^'P P) + ^:^^e/^-"/^5^(J/V.),TV(5,y^P), (7.30) 

where V{P) is the vector (pseudoscalar) nonet matrix. We have g^KK = 4.8 and gj/^K*K = 2.6 x 
10^'^ GeV^^ fitted from the data The gj^^r^s coupling is estimated to be ^-^gj/^K^K-: where 
the deviation comes from SU(3) breaking and the contribution of electromagnetic interactions 
and we estimate gx'Kqs — g<t>KK- 
It is straightforward to obtain 

Abs{J/i;K;7]s) = - / ,^„vsoc^ .o„n.so^ (2^) ^ (.PB - Pi - P2) 

F^it,mn J ^(2) 

t-ml 



2 J (2^)32^1 (2^)32^2 ^ " 2e* • P 

x{-i)[£h^hP3,PA]\ - 2\/2igj/^^^^ +_'^2'' 9K*Kr,A[ 

1 f d^pi d?p2 ,^ ,4,4, .A{^^KJ/il:) 



with 



Abs{K/Jip-K) = - I , , {2-Kfb^(vB-p\-P2) 

^ ' ^' ^ 2 i (27r)32£;i (27r)32S2 2e^ • P 

x(-i)[£*,e:,p3,P4]|2^/2ig^^x ^'^^'T^ /j/^i^*i^4'^j> (7-31) 

L t — J 

^(B ^ J/^K) = ^KfeK:a2/j/^mj/^Ff ^(m2/^)(2e}/^ • ps). (7.32) 

We obtain 

Ai^{J/^ljK; ris) = 2.4 x 10"" GeV, A^^{J/%bK; K) = -6.9 x 10"" GeV, (7.33) 

where use of 02 = 0.28, /j/^ = 405 MeV, ~ 800 MeV and F{t,mK,r}s) — 1 has been made. 
Although the cancellation is incomplete as expected, the FSI contributions (even without the form 
factor suppression) from the B —>■ J/ tpK decay are too small to explain the data (see Tables IVIII 
and IVIlH) . The smallness of these contributions is due to the smallness of the OZI suppressed 
Sj/ippv couphng. 

One possibility of having a large f± is to circumvent the cancellation in S — > VP — > VV 
and B PV — > VV rescatterings. For example, we may consider B — > SA — > cpK* contributions, 
where S and A denotes scalar and axial vector mesons, respectively. We need SA instead of VS, PA 
intermediate states to match the {4)K*)p-uiave quantum numbers. The B — > SA and B AS 
amplitudes may not be similar (note that f^)*^ <^ f^g) and the charge conjugation properties of 
A{^Pi) and A{^Pi) are different. 

To show this, we consider the FSI contributions to f±{VV) from the even-parity charmed 
meson (denoted by D**) intermediate states. The Lagrangian {Cd**d**v) for the interactions of 
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TABLE XI: Short-distance {SD) and long-distance (LD) contributions from intermediate 
states to the decay ampUtudes (in unit of 10"^ GeV) and the polarization fractions for B — > (f)K* . 
LD contributions are calculated with A/j* = mj^* + t/Aqcd and rj = 1.5. 







longitudinal (L) 


parallel ( ) 


perpendicular (±) 






0.63 - 1.98i 


-0.18 + 0.55i 


0.15-0.46i 


j^SD+LD 


-1.65 - 2.07i 


0.50 + 0.58i 


1.29 - 0.42? 






-0.33-O.Oli 


0.00 + 0.00? 







Dsi) 


0.03 + O.OOi 


-1.99 -0.07i 





A{DsiDl 


Dto) 








0.07 + O.OOi 


A{DsiDl 


Dsi) 








1.15 + 0.04ii 




Dso) 








0.02 + O.OOi 




Dsi) 








0.26 + O.OH 


A{DsiDi 




-0.01 - O.OOi 


0.83 + 0.03i 


0.14 + 0.01ii 


A{DsiDi 


Dsi) 


-1.97- 0.07i 


1.83 + 0.07? 


-0.50 -0.02i 


fSD 




0.88 


0.07 


0.05 


jSD+LD 


0.74 


0.06 


0.19 


jcxpt 


0.52 ±0.05 


0.23 ±0.06 


0.25 ±0.04 



the charmed meson Jj" = (0^, l^)i/2 multiplets and vector mesons is similar to 1)3. 15(1 for the s- 
wave charm meson case. In chiral and heavy quark limits, the effective Lagrangian can be expressed 
compactly in terms of superfields [6o| 

Cr)^,D»V = i(3l{Sbv''{V^ - p^,)ba'Sa) + iMSb<y'"'F^,{p),,aSa). (7.34) 

with the superfield 5 = ^^(-0^7/^75 — Dq), which is similar to the superfield H except for the 75 
factor and the phase in front of Dq. Hence, the Lagrangian in terms of the component fields can be 
obtained from Eq. 1)3. 15(1 with qddv, fo^DV-, fD*D*v, gD*D*v, P, ^, L'*-field and D-field replaced 
by go'D'V, fDiD*v, foiDiV, QDiDiV, Pi, Ai, Di-field and -iD^-field, respectively, and with an 
overall negative sign on C. 

The absorptive part contributions of ii? ^ _^ ^j^* j-)*^ ^^^^^ exchanges are given 

by 

Abs iD:QD*o;D*Q or Dsi) = Abs {DsD; Ds or D*) with D*^^^ ^ D(^s)i and ^ D*^^^^, 
Abs {DsiD*o;D*Q or Dsi) = -iAbs (D^D; Ds or D*) with D^^ ^ D^s)i and D(^s) ^ ^(,)o, 
Abs {D:oDr,D*o or Dsi) = iAbs {DsD*-Ds or D^) with Z)^^) ^ D(^s)i and D^s) ^ ^(,)o, 
Abs iDsiDi;D*o or Dsi) = Abs {D*D*-Ds or D*) with D*^^^ ^ D^s)i and D^s) ^ ^(,)o> 



(7.35) 



and the B D**D** factorization amplitudes are 

A(B^D*M = i^V,bV:,ajDt,{rnl-rnl.)F^''Hrnhj 
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A(B- 
A(B - 
A(B 



D*D*) 



.Gf 



-.BD* 



-^VcbVcsaifDsirriDsiFi (m^^ J(2eJ)^^ ■ pb) 



Gf 



2Vr 



BDi 



2A 



{ruB - mD") 



[ruB - mo 



2^fii^al3PBPDi 



(7.36) 



with decay constants and form factors given in 63] • Note that the —i and i factors appearing in 
(|7.35|) originate from the phases differences of -D^^^q and Df^g-^ in the superfields S and H, respectively. 
From these equations, it is now evident that the contributions from B — > DsiDq (pK and 
B D*qDi — > (pK* rescatterings add up instead of cancellation as noticed before. To estimate 
the numerical significance of these contributions, we shall use /3i = 1(~ /?), Ai = lGeV^^(~ A) 
and A£)*^ D^-^ = m£)*^^£,_^-^ + r/AqcD with -q = 1.5 and the monopole form for the form factors 
F{t,m£)*) and F{t,m£)^i). Numerical results are shown in Table IXll The branching ratio is 
j^SD+LD^j^ _^ (f)K*) = 8.4 X 10~^. It is interesting to notice that most of the FSI contributions 
from the intermediate D**D** states add up in A± giving f± ~ 0.2, which is close to data. Note 
that we only consider the contribution from -0^1/2-^1/2 intermediate states. A more detailed study 
including other multiplets such as along this line is worth pursuing. 



F. Comparison with other works 

Alth oug h the analysis of the (f)K* polarization anomaly by Colangelo, De Fazio and Pham 
(CDPHTg] is very similar to ours, their results are quite different form ours in some respects (see 
also j85| for a different analysis based on the Regge theory). First, in our case, the interference 
between short-distance and long-distance contributions is mild because the absorptive long-distance 
contribution is essentially "orthogonal" to the short-distance one (see Table HX)) . By fixing A = 2.3 
GeV, to be compared with our choice of Ad^^d* = mD^^D* + r^AqcD with rj = 0.8, CDP introduced 
a parameter r so that A = Asd + f-^LD with r = corresponding to the absence of rescattering. 



It appears from Fig. 3 of [7g] that there is a drastic SD and LD interference behavior in the CDP 
calculation, namely, the resulting branching ratio is not symmetric with respect to the reflection of r: 
r — > — r. It is not clear to us why the absorptive LD contribution obtained by CDP is not orthogonal 
to the SD one at all. Second, instead of employing a small r factor of order 0.08 to reduce the 
FSI contribution, we use the form factors with different momentum dependence for the exchanged 
particles Dg and D*. As a result, fi ~ 1/2 can be accommodated in our work with a reasonable 
cutoff. Third, we have shown large cancellations occurring in the processes B D*D <j)K* 
and B ^ DsD* ^ (t>K and this can be understood as a consequence of GP and SU(3) [CPS] 
symmetry. We are not able to check the aforementioned cancellation in the CDP work as no details 
are offered to the individual LD amplitudes. Anyway, f± is very small (of order 3%) in our case, 
while CDP obtained f± ~ 0.15. Fourth, CDP adapt a different phase convention for pseudoscalar 
mesons. Considerable efforts have been made in keeping the consistency of phase conventions and 
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reminders on phase subtleties (see, for example, footnote 2 on p. 11) are given in this work. On 
the other hand, as far as the analytic expressions are concerned we do not find any inconsistency 
in the phase convention chosen in the CDP work. Fifth, we apply the U{3) symmetry for the 4> 
meson, while CDP seem to consider only the cjg component of it. Since in both approaches (j) always 
couples to a oi*^ oi*^ pair (see Fig. |^in this work and Fig. 2 in 78|), only the ss component of (f) is 
relevant. Hence, the FSI amplitudes in CDP should be identical to ours up to an overall factor of 
the Clebsch-Gordan coefficient, — 2/-v/6- We do not expect any qualitative difference arising from 
the different treatments of the (p wave function. In particular, the cancellation argument presented 
in Sec. VII. C should be also applicable to their case. 

Working in the context of QCD factorization, Kagan 2a] has recently argued that the lower value 
of the longitudinal polarization fraction can be accommodated in the Standard Model provided that 
the penguin-induced annihilation contributions are taken into account. In general, the annihilation 
amplitudes induced by the {V — A){V — A) operators are subject to helicity suppression. Indeed, 
they are formally power suppressed by order {A.qcd / n^'b)^ (131 • However, the {S — P){S + P) 
penguin operator induced annihilation is no longer subject to helicity suppression and could be 
substantially enhanced. In QCD factorization, the penguin-induced annihilation amplitude contains 
some infrared divergent integral, say X^, arising from the logarithmic endpoint divergence, and 
the logarithmic divergence squared X\. Modelling the logarithmic divergence by 

XA = {l + PAe"^^)ln'^, (7.37) 

with A/j being a hadron scale of order 500 MeV and pa ^ ^, Kagan has shown that the (pK* 
polarization measurements can be accommodated by fitting to the annihilation parameter pA (see 
Fig. 4 of [3]). Note that the penguin- induced annihilation terms being formally of order 1/mf do 
respect the scaling law stated in Eq. (|7.1() . but they can be 0(1) numerically. 

Assuming Kagan's reasoning [22| is correct, it means that although QCD factorization does not 
predict a lower value of the (pK* longitudinal polarization, one may be able to find a value for pA 
within an acceptable range to be able to accommodate the data. 

It is worth mentioning also that the annihilation amplitudes in VV modes are calculable in 
the so-called perturbative QCD approach for hadronic B decays where the end-point divergence is 
regulated by the parton's transverse momentum. In particular, the B — > (pK* decay is studied in 
0, giving /l : /y : f± = 0.75 : 0.13 : 0.11 and B{B- (pK*-) = 16.0 x 10"^. Although fi is 
reduced after the inclusion of nonfactorizable and annihilation contributions, /y _|_ are only of order 
10%. The annihilation contribution indeed helps but it is not sufficient. 

A common feature of the final-state rescattering or a large penguin-induced annihilation contri- 
bution as a mechanism for explaining the observed (pK* anomaly is that the same mechanism will 
also lead to a large transverse polarization in the pK* modes, which is borne out by experiment 
in the decay B~ p~K*^ but not in B^ — > p^K*~ (see Table IVlH) . This has to be clarified 
exp er iment ally. 

An alternative suggestion for the solution of the (pK* anomaly was advocated in js^ that a 
energetic transverse gluon from the b ^ sg chromodipole operator keeps most of its quantum 
numbers except color when it somehow penetrates through the B meson surface and descends 
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to a transversely polarized (j) meson. Sizable transverse components of the B — > tjjK* decay can 
be accommodated by having /y > /_l (see also Table IVnH) . Since the gluon is a flavor singlet, 
this mechanism can distinguish cj) from p, hence it affects B — > cjiK* ^ujK* but not B — > K* p. The 
B —I- ojK* is predicted to be of order 4x 10~^ with large transverse components. On the other hand, 
a recent estimation based on PQCD [s^ seems to indicate the above-mentioned contribution is too 
small to explain the (pK* anomaly. Moreover, the predicted absence of the transverse polarization 
m. B ^ pK* is not consistent with experiment at least for p~ K and p^K*^ modes. 

VIII. DISCUSSION AND CONCLUSION 

In this work we have studied the effects of final-state interactions on the hadronic B decay rates 
and their impact on direct CP violation. Such effects are modelled as soft final-state rescattering 
of some leading intermediate two-body channels and can be classified according to the topological 
quark diagrams. It amounts to considering the one-particle-exchange processes at the hadron level 
for long-distance rescattering effects. Our main results are as follows: 

1. The color-suppressed neutral modes such as B^ D , K^t:^ can be substan- 
tially enhanced by long-distance rescattering effects, whereas the color-allowed modes are 
not significantly affected by FSIs. 

2. All measured color-suppressed charmful decays of B^ into D^*^^tt^ , D^r], D'^uj and D^p^ are 
significantly larger than theoretical expectations based on naive factorization. The rescat- 
tering from B — > {Dtt, D*p} Dtt contributes to the color-suppressed VF-emission and 
ly-exchange topologies and accounts for the observed enhancement of the D^-k^ mode with- 
out arbitrarily assigning the ratio of 02/01 a large magnitude and strong phase as done in 
many previous works. 

3. The branching ratios of all penguin-dominated B irK decays are enhanced via final-state 
rescattering effects by (30-40)%. CP asymmetry in the K~7r~^ mode is predicted by the 
short-distance approach to be ~ 4%. It is modified by final-state rescattering to the level of 
14% with a sign flip, in good agreement with the world average of —0.11 it 0.02. Given the 
establishment of direct CP violation in B^ — > K~it~^, a model-independent relation assuming 
SU(3) implies that the vr+vr^ mode has direct CP asymmetry of order 40% with a positive 
sign. An isospin sum-rule relation involving the branching ratios and CP asymmetries of the 
four B Kit modes indicates the presence of electroweak contributions and/or some New 
Physics. Better measurements are needed to put this on a firmer footing. 

4. If only the absorptive part of final-state rescattering is considered, it cannot explain the 
observed enhancement of the rate, whereas the predicted vr^vr will be too large by 
a factor of 2 compared to experiment. The dispersive part of rescattering from DD — > vrvr 
and TTTT — > TTTT via meson annihilation (i.e. the so-called vertical VF-loop diagram V) which 
interferes destructively with the short-distance amplitude of B^ vt+tt" can reduce vr+vr" 
and enhance substantially. 
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5. For tree-dominated B vrvr decays, it is known that the predicted direct CP asymmetry 
in the tt^tt^ mode by the short-distance approach is small with a negative sign. We have 
shown that its sign is flipped by final-state rescattering and its magnitude is substantially 
enhanced. Direct CP violation in i? — > vr'^vr'' is predicted to have a sign opposite to that of 
vr^vr", in contrast to the predictions based on perturbative QCD (PQCD). Hence, even a 
sign measurement of direct CP asymmetry in tt^tt^ can be used to discriminate between the 
FSI and PQCD approaches for CP violation. 

6. CP partial rate asymmetry in — > vr^vr'^ is very small in the SM and remains so even 
after the inclusion of FSIs, with a magnitude less than 1 percent. Since CP asymmetry is 
caused by the electroweak penguin but not the QCD penguin (up to isospin violation), this 
is a good mode for searching New Physics. 

7. For tree-dominated B ^ pir decays, we showed that (i) the color-suppressed p^Tr^ mode is 
slightly enhanced by rescattering effects to the order of 1.3x 10~^, which is consistent with the 
weighted average of the experimental values. However, the discrepancy between BaBar and 
Belle for this mode should be clarified soon. It should be stressed that direct CP violation in 
this mode is significantly enhanced by FSI from around 1% to 60%. (ii) Direct CP violation 
in the p~^tt~ mode is greatly enhanced by FSI from the naive expectation of ~ —0.01 from 
the short-distance approach to the level of —0.42, in agreement with BaBar and Belle. 

8. As for the intriguing B (pK* polarization anomaly, the longitudinal polarization fraction 
can be significantly reduced by the rescattering contribution from the intermediate DDs 
states. However, no sizable perpendicular polarization is found owing mainly to the large 
cancellations among various contributions from intermediate DgD* and D*D states. Conse- 
quently, our result for the perpendicular polarization fraction is different from a recent similar 



analysis in [7g] . Final-state rescattering effects from this particular set of states seem not be 



able to fully account for the polarization anomaly observed in i? — > (j)K* and FSI from other 



intermediate states and/or some other mechanism e.g. the penguin-induced annihilation |27l |. 
may have to be invoked. Given the fact that both Belle and BaBar observe large phases in 
various polarization amplitudes (see Table IVTH) . FSI may still provide a plausible explana- 
tion. In any case, our conclusion is that the small value of the longitudinal polarization in 
VV modes cannot be regarded as a clean signal for New Physics. 

Needless to say, the calculation of final-state rescattering effects in hadronic B decays is rather 
complicated and very much involved and hence it suffers from several possible theoretical uncer- 
tainties. Though most of them have been discussed before, it is useful to make a summary here: 

(i) The strong form factor F[t) for the off-shell effects of the exchanged particle. It is param- 
eterized as in Eq. (|3.22jl by introducing a cutoff scale A. Moreover, we write A = m-exc + ??Aqcd 
[cf. Eq. ()3.40|) ]. where the parameter rj is expected to be of order unity and can be determined 
from the measured rates. For a given exchanged particle, the cutoff varies from process to process. 
For example, rir, = 2.1, 0.69, 1.6 for B — > Dtt, Ktt, pir decays, respectively. To see the sensitivity 
on the cutoff, we have allowed 15% error in the QCD scale Aqcd- Another important uncertainty 
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arises from the momentum dependence of the form factor F{t). Normally it is assumed to be of the 
monopole form. However, the analysis of the (pK* decays prefers to a dipole behavior for the form 
factor F{t,mD*) appearing in the VexcVP or VcxcW vertex with the exchanged vector particle 
V^xc and a monopole dependence for other form factors. A more rigorous study of the momentum 
dependence of the strong form factor is needed. It turns out that CP asymmetries and the decay 
rates especially for the color-suppressed modes are in general sensitive to the cutoff scale. This 
constitutes the major theoretical uncertainties for long-distance rescattering effects. 

(ii) Form factors and decay constants. Model predictions for the form factors can vary as much 
as 30% which in turn imply large uncertainties in the branching ratios. In particular, the B VV 
decays are fairly sensitive to the difference between the form factors Ai and A2- For form factor 
transitions in this work we rely on the covariant light-front model which is favored by the current 
data. 

(iii) Strong couplings of heavy mesons and their SU(3) breaking. We have applied chiral and 
heavy quark symmetries to relate various strong couplings of heavy mesons with the light pseu- 
doscalar or vector mesons. It is not clear how important are the chiral and l/mg corrections. For 
SU(3) breaking effects, we have assumed that they are taken into account in the relations given in 
Eqs. KW\ and dZH). 

(iv) The real part of the long-distance contribution which can be obtained from the disper- 
sion relation (|3.39() . Unlike the absorptive part, the dispersive contribution suffers from the large 
uncertainties due to some possible subtractions and the complication from integrations. For this 
reason, we have ignored it so far. However, in order to resolve the discrepancy between theory 
and experiment for tt+tt" and B^ — > vr'^Tr*^, we have argued that it is the dispersive part of 
long-distance rescattering of DD vrvr and vrvr vrvr via meson annihilation that accounts for the 
suppression of the vr^vr^ mode and the enhancement of 

(v) CKM matrix elements and 7. Direct CP asymmetry is proportional to sin 7. All numbers 
in the present work are generated by using 7 = 60°. 

(vi) Intermediate multi-body contributions which have not been considered thus far and may 
have cancellations with the contributions from two-body channels. 

Finally, it is worth remarking that the final-state rescattering effects from intermediate charmed 
mesons as elaborated in the present work has some similarity to the long-distance "charming pen- 
guin" effects advocated in the literature 32 , |3J] . The relevance of this long-distance effect has 



been conjectured to be justified in the so-called soft collinear effective theory [89|. Indeed, we have 
pointed out that some of the charmed meson loop diagrams in the decays e.g. B — > vrvr, ttK, (pK* 
will manifest as the long-distance cc penguins. However, we have also considered FSIs free of 
charming penguins. For example, long-distance rescattering effects in S — > Dtt and in B~ — > ■k~tt^ 
(also p~ p^) decays have nothing to do with the charming penguin effects. In other words, our 
systematic approach for FSIs goes beyond the long-distance charming penguin mechanism. 
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APPENDIX A: USEFUL FORMULA 



Under the integration, the covariant integrals 



(27r) 5 {pB -Pl -P2)f{t) X {pif,, pif^piu, Pl^iPluPla} 



(Al) 



(27r)32Ei (27r)32^2 

can only be expressed by external momenta ^3,^4 with suitable Lorentz and permutation structures. 
Hence we can express these plf^, pi^piu, Pi^PiuPia in terms for P = P3 + p^ and q = p^ — p^ as in 
Eq. ()7.12|) within the integration. By contracting the left-hand side of Eq. 1)7. 12() with P^, q^ and 
9^11/, we are able to solve for these A^^^ and obtain 
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(3) 

where the lower left part of the symmetric inverse matrix for the Aj case is not shown explicitly. 
Note that the t dependence of A^ is in the column matrices in the right-hand side of the above 

(i) 

equations and t always appears in the numerators of A - in a polynomial form. 



APPENDIX B: THEORETICAL INPUT PARAMETERS 



In this Appendix we summarize the input parameters used in the present paper. For decay 
constants we use 

= 132 MeV, fx = 160 MeV, fo = 200 MeV, fo* = 230 MeV, fa, = -205 MeV. (Bl) 

Note that a preliminary CLEG measurement of the semileptonic decay — > n'^i' yields f n+ = 
(202 zb 41 zb 17) MeV |90|. For form factors we follow the covariant light-front approach 63 1. For 
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purpose of comparison, we list some of the form factors at = used in the main text: 

Fo^^(O) = 0.25, Fo^^(O) = 0.35, A^P{0) = 0.28, A^^* (0) = 0.31 (B2) 
to be compared with 

Fo^'^(O) = 0.28 ±0.05, Fo^^(O) = 0.34 ±0.05, ^^''(0) = 0.37 ± 0.06, (0) = 0.45 ± 0.07 

(B3) 

employed in 0| and 

F(f^(0) = 0.33, F(f^(0) = 0.38, A^^(O) = 0.28, v4^^*(0) = 0.32 (B4) 



used by Bauer, Stech and Wirbel |9l|. The most recent light-cone sum rule analysis yields 92 1 

Fo^'^(O) = 0.258 ± 0.031, F^^ (0) = 0.331 + 0.041 + 0.025(5„, (B5) 

with being defined in |92l]. The B tt^tt^ data favors a smaller F^^{0) of order 0.25, while 
B ^ pn measurement prefers to small form factors for B ^ p transition. 

For the quark mixing matrix elements, we use A = 0.801 and A = 0.2265 [t^] in the Wolfenstein 
parametrization of the quark mixing angles j^^. The other two Wolfenstein parameters p and rj 
obey the relations f] = i?sin7 and p = Rcosj, where R = \ VudVub/ i^cd.^cb)\^ P = ~ ^^/'^)P 
f] = {1 — A^/2)r/. We take 7 = 60° and R = 0.39 in this work for calculations. For current quark 
masses, we use mh{mi,) = 4.4 GeV, mc{mb) = 1.3 GeV, ms(2.1 GeV) = 90 MeV and niq/ms = 0.044. 

The physical strong coupling constants are 

(7p.^ = 6.05 ±0.02, 5^,^^ = 4.6, ^D'Dtt = 17.9 ±0.3 ±1.9, 
gv = 5.8, /3 = 0.9, A = 0.56GeV"\ (B6) 

where gv, P and A (not to be confused with the parameter A appearing in the parametrization of 
quark mixing angles) are the parameters in the effective chiral Lagrangian describing the interac- 
tions of heavy mesons with low momentum light vector mesons. 
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